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A 

In  this  dissertation  I  investigate  equatorial  ocean  dynamics 
using  analytical  and  simple  numerical  models.  The  steady  circulation 
Is  first  examined  using  a  linear  equatorial  beta-plane  formulation  In 
which  continuous  density  stratification  Is  a  key  element.  The  basin  Is 
horizontally  unbounded  and  forced  by  either  winds  or  a  surface  tempera¬ 
ture  boundary  condition.  Two  vertically  diffusive  models  are  proposed 
to  examine  the  effects  of  different  turbulence  parameterlzatlons.  One 
Is  characterized  by  temperature  dissipation  proportional  to  temperature 
Itself,  the  other  has  biharmonic  friction  proportional  to  a  fourth 
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derivative  of  velocity.  Both  lead  to  equations  for  the  steady  state 
that  are  In  many  ways  analogous  to  the  well-known  equations  governing 
equatorial  wave  motion.  The  analogies  are  exploited  both  In  obtaining 
analytical  solutions  and  In  describing  the  flow  dynamics. 

Steady  zonal  winds  drive  geos trophic  flows  such  as  the 
Equatorial  Undercurrent  and  South  Equatorial  Current  as  well  as  equa¬ 
torial  upwelllng  and  a  weakening  of  the  equatorial  thermocllne.  Exami¬ 
nation  of  the  dynamics  reveals  that  In  a  stratified  ocean  equatorial 
Ekman  layers  are  well  behaved  due  to  the  presence  of  meridional  baro- 
cllnic  pressure  gradients.  These  pressure  gradients  are  not  geostrophlc- 
ally  balanced  and  therefore  Imply  an  error  of  0(10%)  or  greater  In  the 
geostrophlc  calculation  of  near  surface  zonal  currents  from  the  density 
field.  Furthermore,  the  geostrophlc  balance  Is  modified  near  the  equator 
to  allow  for  direct  wind  forcing  and  vertical  diffusion  of  zonal  momentum. 
Nonlocal  wind  effects  are  Important  In  the  model  results  suggesting  a 
possible  dynamical  connection  between  maximum  westward  winds  in  the 
central  equatorial  Pacific  and  maximum  undercurrent  speeds  in  the  eastern 
Pacific.  Nonlinearity  Is  Important  In  the  zonal  momentum  balance,  but 
Its  effects  on  the  undercurrent  can  be  qualitatively  predicted  by  quad¬ 
ratic  corrections  to  the  linear  model.  The  above  results  apply  to  either 
turbulence  parameterization  scheme;  only  the  magnitude  of  various  flow 
quantities  changes. 

Large  scale  meridional  winds  or  zonal  winds  with  meridional 
shear  excite  a  frictional  response  that  is  too  shallow  to  displace  deeper 
geostrophlc  zonal  flows  off  the  equator.  Thus  either  time  dependence 
or  strong  nonlinearity  must  be  invoked  to  account  for  observations  of  an 
undercurrent  structure  asymmetric  about  the  equator.  Surface  zonal  tern- 


perature  gradients  drive  a  current  system  similar  to  that  driven  by  zonal 
winds  suggesting  that  the  equatorial  mixed  layer  will  have  a  significant 
Impact  on  the  deeper  circulation. 

The  Interaction  of  prescribed  geostrophlc  zonal  mean  currents 
In  an  equatorial  ocean  with  free,  neutrally  stable  Kelvin  and  Inertlo- 
gravlty  waves  Is  Investigated  next.  The  equations  are  formulated  for  a 
two  layer  ocean,  then  solved  either  analytically  or  by  simple  numerical 
methods.  It  Is  found  that  the  zonal  velocity  of  Inertlo-gravlty  waves  Is 
strongly  affected  In  regions  where  mean  currents  have  a  large  meridional 
shear.  This  coupling  can  be  Interpreted  as  mean  current  meanders  simi¬ 
lar  to  those  observed  during  GATE.  Conversely,  sea  level  fluctuations 
proportional  to  wave  pressure  are  not  as  greatly  influenced  by  mean 
currents.  These  results  may  be  of  Importance  when  attempting  to  fit 
equatorial  wave  theories  to  sea  level  observations  on  one  hand  and 
current  measurements  on  the  other. 

Solutions  In  the  Kelvin  wave  range  indicate  that  at  low  frequen¬ 
cies,  the  presence  of  mean  currents  leads  to  a  small  meridional  wave 
velocity  and  a  Doppler  shift  in  frequency.  This  latter  effect  may  in¬ 
fluence  the  speed  at  which  the  equatorial  ocean  adjusts  to  transient 
forcing.  At  higher  frequencies,  the  Kelvin  wave  becomes  more  gravity¬ 
like  and  may  appear  In  current  records  as  varicose  meanders  of  the  mean 
flow. 
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CHAPTER  I 

INTRODUCTION 

The  equatorial  ocean  has  been  a  continuing  source  of  Interest  to 
oceanographers  since  the  first  modern  observations  of  the  Pacific 
Equatorial  Undercurrent  in  1952  (Montgomery  and  Stroup,  1962).  This 
strong  eastward  current  flows  below  a  generally  westward  wind  drift  and 
extends  no  further  than  2°  latitude  from  the  equator.  Similar  under¬ 
currents  were  discovered  shortly  thereafter  in  both  the  Indian  Ocean  dur¬ 
ing  the  Northeast  Monsoon  (Taft  and  Knauss,  1967)  and  the  Atlantic  (Neu¬ 
mann,  1960).  Even  more  subsurface  structure  emerges  from  recent  analyses. 
Taft  et  al.  (1973)  summarize  evidence  for  a  persistent  westward  flow  be¬ 
low  the  permanent  undercurrents;  Tsuchlya  (1975)  and  Cochrane  et  al. 
(1979)  describe  eastward  subsurface  countercurrents  flanking  this  west¬ 
ward  flow,  luyten  and  Swallow  (1976)  and  Luyten  and  Eriksen  (1980)  have 
observed  zonal  jets  trapped  near  the  equator  down  to  4000  m.  In  addition 
to  this  rich  spatial  structure,  there  is  a  full  spectrum  of  geophysically 
interesting  temporal  fluctuations.  For  example,  evidence  of  a  sea  level 
signal  at  periods  of  4-5  days  exists  in  the  equatorial  Pacific  (Wunsch 
and  Gill,  1975).  Undercurrents  have  meandered  about  the  equator  at 
periods  of  0(10  days)  In  both  the  Atlantic  (DUIng  et  al.,1975)  and 
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Pacific  (Titov,  1977).  At  longer  time  scales,  the  equatorial  Indian  and 
Pacific  Oceans  respond  to  seasonal  monsoon  forcing.  Interannual ly,  large 
scale  coupled  ocean-atmosphere  disturbances  known  as  El  NiPio  occur. 

Daunted  by  the  variety  and  complexity  of  oceanic  phenomena, 
Melville  wrote  In  1851,  "...the  secrets  of  the  currents  in  the  seas  have 
never  yet  been  divulged,  even  to  the  most  erudite  research...".  But 
through  technological  advances,  systematic  observations  and  the  study  of 
dynamical  models,  we  have  in  the  last  30  years  gained  a  wealth  of  Insight 
Into  the  peculiar  physics  of  the  equatorial  ocean.  For  instance.  It  has 
been  possible  to  Identify  some  of  the  important  mechanisms  that  generate 
and  maintain  the  undercurrent  system.  Likewise,  linear  equatorial  wave 
theory  (Moore  and  Philander,  1977)  Is  a  unified,  elegant  paradigm  for 
describing  a  wide  range  of  time  dependent  phenomena.  From  this  body  of 
research,  two  fundamental  concepts  emerge.  First,  the  equatorial  ocean 
Is  unique  dynamically  because  horizontal  Coriolis  forces,  of  primary  im¬ 
portance  in  mid-latitude  flows,  vanish  at  the  equator.  Second,  vertical 
structure  near  the  equator  depends  intimately  on  density  stratification. 

From  a  broader  perspective,  progress  In  equatorial  oceanography 
moves  us  closer  to  understanding  the  ocean-atmosphere-biosphere  complex 
and  Its  Impact  on  human  affairs.  For  example,  in  sqme  years  a  sudden 
warming  of  waters  off  the  Peruvian  coast  devastates  the  anchovy  fishery 
and  produces  torrential  rains  in  a  normally  arid  region.  These  events 
known  as  El  Nlfio  appear  to  be  related  to  large  scale  wind  fluctuations 
thousands  of  kilometers  to  the  west  In  the  Central  Pacific  (McCreary, 
1977).  A  better  understanding  of  such  events  could  alleviate  many  of 
their  adverse  social,  political  and  economic  effects.  Also  primary 
productivity,  the  base  of  the  food  chain.  Is  strongly  linked  to 


3 


upwelling  which  brings  nutrient  rich  water  to  the  surface.  Nowhere  out¬ 
side  of  coastal  regions  are  there  as  Intense  large  scale  vertical  veloci¬ 
ties  as  at  the  equator.  Knowledge  of  the  factors  that  control  the  spatial 
and  temporal  aspects  of  equatorial  upwelling  (and  other  vertical  trans¬ 
port  processes  e.g.,  mixing)  Is  thus  of  paramount  Importance  In  under¬ 
standing  the  productivity  of  this  region.  Finally  In  the  tropics.  Inci¬ 
dent  solar  energy  exceeds  outgoing  radiative  flux.  The  excess  is  trans¬ 
ported  poleward  by  atmospheric  and  oceanic  processes  that  create  weather 
and  climatic  variability.  Recent  studies  (e.g.  Shukla,  1975;  Reiter, 

1978)  suggest  that  equatorial  current  and  temperature  structures  are  key 
elements  In  the  generation  of  this  variability  not  only  In  the  tropics  but 
also  at  temperate  latitudes. 

In  this  dissertation  I  examine  the  physics  of  steady  state  and 
time  variable  flows  In  linearized,  stratified,  equatorial  beta-plane 
models.  Chapter  II  begins  with  a  discussion  of  observations  pertinent  to 
the  steady  models  and  a  brief  review  of  previous  theories.  I  then  develop 
two  models  for  steady  flow  In  a  continuously  stratified  ocean  to  illustrate 
the  boundary  layer  nature  of  the  equatorial  regime  and  the  transition 
from  equatorial  to  mid-latitude  dynamics.  The  basic  conclusions  are  that 
(a)  stratification  removes  singularities  In  the  equatorial  extension  of 
Ekman  layers,  (b)  mid-latitude  geostrophy  Is  modified  near  the  equator  to 
Include  vertical  diffusion  of  zonal  momentum  and  (c)  zonal  winds  near  the 
equator  enter  the  ocean  not  through  the  Ekman  layer  but  through  direct 
projection  onto  geostrophlc  currents.  It  Is  also  shown  that  (d)  stratifi¬ 
cation  Is  crucial  In  setting  the  scales  of  various  flow  parameters  and  that 
(e)  barocllnlc  zonal  pressure  gradients  are  essential  for  the  maintenance 
of  vertical  structure.  Moreover  these  results  apply  for  different 
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turbulence  parameterizations.  Chapter  III  Is  a  study  of  the  Inter¬ 
action  of  free  waves  and  prescribed  zonal  shear  flows  like  the  under¬ 
current  and  westward  South  Equatorial  Current.  Stratification  in  this 
model  Is  approximated  by  a  single  jump  In  density  near  the  surface.  The 
important  results  of  this  section  are  (a)  Internal  equatorial  waves  can 
cause  mean  current  meanders  like  those  observed  In  the  Pacific  and 
Atlantic  and  (b)  mean  flow,  depending  on  its  direction,  can  either  en¬ 
hance  or  retard  the  oceanic  response  time  to  transient  wind  forcing  near 
the  equator. 


CHAPTER  II 

STEADY  STATE  MODELS 
A.  Observational  Background 

Figure  1  presents  the  structure  of  zonal  flow  and  temperature  on 
a  meridional  plane  in  the  western  Pacific  from  Hisard  et  al.  (1970).  The 
basic  pattern  is  representative  of  steady  conditions  at  other  longitudes 
in  the  Pacific  and  Atlantic  Oceans.  It  consists  of  (a)  an  eastward 
Equatorial  Undercurrent  (EUC)  in  the  upper  thermocline  at  speeds  of 
0(100  cm  sec”1),  (b)  a  westward  surface  flow,  the  South  Equatorial 
Current  (SEC)  at  speeds  typically  of  0(50  cm  sec-1),  (c)  a  weaker  west¬ 
ward  Equatorial  Intermediate  Current  (EIC)  below  the  EUC  and  (d)  subsur¬ 
face  countercurrents  (SSCC)  flanking  the  EIC  and  at  this  longitude  con¬ 
tiguous  to  the  shallower  EUC.  The  North  Equatorial  Countercurrent  (NECC) 
Is  visible  at  the  northern  extreme  of  the  section.  All  these  zonal  cur¬ 
rents  are  In  geostrophlc  balance  with  the  density  field  which,  to  zeroth 
order,  Is  determined  by  the  temperature  field  of  Figure  lb.  In  particu¬ 
lar,  the  EUC  Is  associated  with  a  weakening  of  the  thermocline  near  the 
equator  and  the  SSCC's  are  in  balance  with  the  plunging  10°  isotherm 
located  symmetrically  at  38N  and  3°S.  A  secondary  eastward  flow  which 
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Figure  1.  (a)  Meridional  section  of  zonal  velocity  at  170°E.  Labelled 

currents  are  South  Equatorial  Current  (SEC),  Equatorial  Undercurrent 
(EUC),  Equatorial  Intermediate  Current  (EIC),  North  Equatorial  Counter- 
current  (NECC)  and  subsurface  countercurrents  (SSCC).  The  current 
labelled  AG  Is  an  ageostrophlc  flow.  Regions  of  westward  flow  are 
stippled,  (b)  Corresponding  meridional  section  of  temperature,  (after 
Hlsard,  et  al.,  1970). 
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lies  between  the  SEC  and  EUC  Is  peculiar  to  those  regions  which  have  a 
deep  mixed  layer.  It  is  neither  permanent  nor  in  geostrophlc  balance 
(Hlsard  et  al.  ,  1970). 

The  source  of  momentum  for  the  EUC  Is  an  eastward  pressure 
gradient  force  set  up  by  the  Trade  Winds.  With  only  weak  deflecting 
forces  near  the  equator,  water  simply  flows  eastward  down  this  pressure 
gradient  below  the  action  of  the  wind  stress.  Using  hydrographic  data, 
Knauss  (1966)  In  the  Pacific  and  Katz  et  al.  (1977)  in  the  Atlantic  have 
measured  this  force  to  be  about  5  x  10-s  dyne  cm-3  at  the  surface  falling 
to  half  that  value  at  the  depth  of  the  EUC  core.  The  signal  is  eventu¬ 
ally  swamped  by  noise  at  greater  depths. 

Cromwell  (1953)  and  later  Fofonoff  and  Montgomery  (1955)  suggest¬ 
ed  a  meridional  circulation  associated  with  flow  patterns  like  those  of 
Figure  1.  Upwelllng,  as  inferred  from  doming  isotherms  and  nutrient 
concentrations  near  the  equator,  feeds  a  surface  meridional  Ekman  diver¬ 
gence.  This  upwelling  and  a  downwelling  inferred  from  downward  trending 
isotherms  below  the  core  of  the  EUC  requires  meridional  convergence  of 
mass  in  the  thermocllne.  The  convergence  is  geostrophically  balanced 
by  the  same  pressure  gradient  that  drives  the  EUC.  Knauss  (1966)  has 
measured  meridional  velocities  at  150°W  of  0(1-10  cm  sec-1)  and  by 
scaling  arguments  deduces  vertical  velocities  of  0(1-5  x  10-3  cm  sec-1). 

The  Indian  Ocean  is  dominated  by  seasonal  monsoons  and  mean  zonal 


westerly  winds  rather  than  steady  easterlies.  Consequently  only  during 
the  Northeast  Monsoon  will  the  circulation  resemble  the  persistent 
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Pacific  and  Atlantic  patterns.  More  comprehensive  surveys  of  these 
flow  features  in  all  equatorial  oceans  can  be  found  in  Knauss  (1963); 
Philander  (1973b);  Leetmaa,  McCreary  and  Moore  (1980);  Tsuchiya  (1975); 
Cochrane  et  al.  (1979)  and  Taft  et  al.  (1974). 

The  recent  observations  of  Luyten  and  Swallow  (1976)  in  the 
Indian  Ocean  and  Luyten  and  Eriksen  (1980)  in  the  western  Pacific  add  a 
new  dimension  to  investigations  of  the  equatorial  circulation.  They 
measured  several  zonal  currents  trapped  to  the  equator  down  to  4000 
meters.  The  data  fail  to  resolve  the  time  scales  of  these  stacked  jets, 
but  they  must  be  at  least  several  months.  Furthermore  it  is  not  clear 
how  or  if  they  are  related  to  the  winds  or  shallower  currents.  It  is 
evident  though  that  the  short  vertical  wavelengths  of  0(100' s  of  meters) 
imply  a  much  higher  degree  of  baroclinicity  than  at  mid-latitudes. 

B.  Previous  Theories 

In  this  section  I  will  discuss  the  development  of  key  dynamical 
concepts  deduced  from  previous  theories.  I  will  not  attempt  a  comprehen¬ 
sive  review  since  there  are  several  published  surveys  of  equatorial 
models,  notably  Philander  (1973b),  Gill  (1975),  Moore  and  Philander  (1977) 
and  Leetmaa,  McCreary  and  Moore  (1980).  I  will  furthermore  restrict  dis¬ 
cussion  to  those  models  in  the  oceanic  parameter  range  rather  than  those 
more  appropriate  to  laboratory  studies,  though  the  latter  can  contribute 
to  our  understanding  of  oceanic  processes. 

The  primary  intent  of  nearly  all  steady  state  equatorial  models 
has  been  to  provide  a  dynamical  description  of  the  EUC  since  this  was  the 
first  subsurface  flow  discovered  and  since  It  is  typically  the  strongest 
of  all  equatorial  currents.  One  feature  common  to  all  these  models  is 
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that  they  are  formulated  on  the  equatorial  beta-plane  i.e.  the  Coriolis 
parameter  is  everywhere  equal  to  a  constant  times  latitude.  This  implies 
that  the  equatorial  zone  is  a  singular  region  of  the  global  circulation 
because  rotational  forces  cannot  dominate  the  flow  there.  Robinson 
(1960)  points  out  that  close  to  the  equator,  a  number  of  processes  can 
compensate  for  the  diminished  Coriolis  forces.  A  model  that  includes  at 
least  some  of  these  may  give  a  reasonably  good  approximation  of  the  flow 
patterns  and  scales.  Accordingly  a  great  many  models  which  cover  a  wide 
range  of  parameter  space  have  been  proposed.  It  is  instructive  to  dis¬ 
tinguish  between  those  that  assume  constant  density  and  those  that  include 
continuous  vertical  density  stratification. 

1.  Constant  Density  Models 

Because  the  EUC  is  a  near  surface  feature,  the  depth  of  the 
ocean  in  constant  density  models  is  taken  to  be  the  depth  of  a  deep  mixed 
layer  like  that  in  the  western  Pacific.  The  assumption  is  that  a  sharp 
thermocline  at  the  base  of  the  mixed  layer  completely  inhibits  communica¬ 
tion  with  the  deeper  ocean.  This  is  not  entirely  justifiable  (e.g.  Knox, 
1976)  but  it  typically  leads  to  a  set  of  mathematically  tractable  equations. 
The  models  of  Stommel  (1960)  and  Charney  (1960)  exemplify  and  Illustrate 
the  limitations  of  this  type  of  formulation.  In  these,  the  depth  Inte¬ 
grated  pressure  gradient  is  in  balance  with  an  easterly  wind  stress.  Via 
the  hydrostatic  relation  this  approximates  the  observed  rise  in  sea  level 
to  the  west  in  both  the  Atlantic  and  Pacific  Oceans.  Other  zonal  varia¬ 
tions  are  assumed  small  compared  to  meridional  variations  so  that  flow 
can  be  examined  solely  on  the  meridional  plane. 

In  Stommel 's  linear  calculation  the  pressure  gradient  is  able  to 
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drive  an  eastward  flow  below  a  westward  wind  drift  near  the  equator 
(Figure  2a).  However  the  scales  are  critically  dependent  on  the  magni¬ 
tude  of  a  vertical  eddy  viscosity.  In  fact  it  is  not  possible  to  adjust 
this  parameter  to  give  realistic  width  and  velocity  scales  simultaneously. 
Furthermore,  Stommel  chose  a  no  stress  rather  than  no  slip  condition  at 
the  bottom  of  his  model  ocean.  Chamey  used  the  latter  condition  in  a 
linear  calculation  and  found  no  eastward  flow  near  the  equator  in  the 
presence  of  easterly  winds. 

The  EUC  in  Chamey1  s  model  depends  on  nonlinearity  (Figure  2b). 
Westward  momentum  fed  in  through  the  surface  by  the  winds  is  carried 
poleward  resulting  in  a  net  gain  of  eastward  momentum  near  the  equator. 
Such  a  tendency  is  predicted  by  the  calculation  of  quadratic  corrections 
to  the  linear  model  (Robinson,  1966).  This  gain  of  eastward  momentum  is 
also  consistent  with  the  tendency  to  conserve  absolute  vorticity  near 
the  equator  in  the  presence  of  an  eastward  pressure  force  (Fofonoff  and 
Montgomery,  1955;  Cane,  1979).  A  particle  as  it  approaches  the  equator 
loses  planetary  vorticity  which  must  be  compensated  for  by  a  gain  of 
relative  vorticity.  The  sign  of  this  gain  is  such  that  the  current 
accelerates  eastward.  However  if  the  friction  coefficient  is  too  small 
in  Charney's  model  ( <14  cm2  sec"1),  solutions  become  unstable. 

The  meridional  circulation  in  these  models  is  similar  to  that 
deduced  from  observations.  In  particular,  for  westward  winds  there  is 
upwelllng  near  the  equator  and  surface  divergence  while  at  depth  there 
is  a  convergence  of  mass.  The  distribution  of  temperature  is  not  pre¬ 
dicted  since  the  entire  layer  is  of  a  single  density. 

While  these  models  shed  light  onto  the  dynamics  of  an  equatorial 
fluid,  they  have  some  serious  drawbacks.  Already  mentioned  are  their 
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Figure  2.  Profiles  of  zonal  velocity  at  the  equator  computed 
from  the  constant  density  models  of  (a)  Stommel  (1960)  and 
(b)  Charney  (1960).  Decreasing  values  of  c  In  the  Charney 
calculation  Imply  stronger  nonlinearity.  Negative  veloci¬ 
ties  are  to  the  west. 
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sensitivity  to  the  size  of  eddy  coefficients  and  type  of  bottom  boundary 
conditions  as  well  as  an  Inability  to  predict  the  temperature  field.  The 
lack  of  longitudinal  variation  Is  not  reasonable.  For  example,  zonal 
variability  is  present  In  the  winds  (Wyrtki  and  Meyers,  1975;  Hellerman, 
1967),  surface  temperature  (Colin  et  al.,  1971 )t  current  speeds,  transports 
and  depth  of  the  EUC  (Knauss,  1960).  Knauss  (1966)  computes  significant 
downstream  advectlon  of  zonal  momentum  In  the  EUC.  Moreover  these  models 
tell  us  nothing  of  the  deeper  EIC  and  SSCC.  Their  crucial  limitation 
though  is  that  they  place  the  EUC  in  the  mixed  layer,  not  in  the  thermo- 
cllne  where  It  is  observed.  Shear  profiles  in  a  constant  density  ocean 
cannot  be  In  geostrophlc  balance  whereas  measurements  show  that  all  the 
permanent  zonal  equatorial  flows  are.  As  such  these  models  are  at  best 
applicable  to  ageostrophlc  flow  like  that  observed  at  the  base  of  the 
deep  mixed  layer  in  the  western  Pacific  (Figure  1).  To  understand  flow 
in  the  thermocline,  clearly  we  require  a  theory  that  includes  continuous 
stratification. 

2.  Stratified  Models 

Even  before  the  discovery  of  the  EUC,  Sverdrup  (1947)  proposed  a 
theory  of  the  steady  circulation  in  a  stratified  equatorial  ocean.  He 
derived  relationships  between  various  depth- Integrated  flow  quantities 
and  surface  wind  stress.  While  such  a  theory  can  provide  no  details  of 
barocllnlc  structure,  Arthur  (1960)  showed  that  Sverdrup's  calculations 
were  consistent  with  the  notion  of  an  eastward  subsurface  flow  Imbedded 
In  a  westward  surface  wind  drift.  Yoshida  (1959)  was  the  first  to  ex¬ 
plicitly  model  the  EUC  In  the  thermocline.  Working  with  a  linearized 
system  of  equations,  he  deduced  a  geostrophlc  balance  for  the  zonal  flows 
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as  well  as  the  dependence  of  certain  scales,  e.g.  width  and  velocity,  on 
stratification.  However  his  results  require  a  large  number  of  a  priori 
assumptions  for  computational  expedience.  For  Instance,  he  neglects  the 
vertical  diffusion  of  meridional  momentum  thus  filtering  out  Ekman  layers. 
Veronls  (I960)  attempted  a  nonlinear  calculation  on  a  zonal  plane  along 
the  equator  by  using  perturbation  expansions  In  latitude.  Though  In 
some  respects  temperature  and  zonal  velocity  resembled  observed  fields, 
the  system  of  equations  could  not  be  closed.  The  perturbation  expansions 
generated  more  unknowns  than  physical  constraints. 

Pedlosky  (1969)  discussed  axl symmetric  flow  of  a  linear  strati¬ 
fied  fluid  between  two  concentric  spheres  in  a  parameter  range  appropri¬ 
ate  to  the  ocean.  For  an  Imposed  surface  heat  flux  and  velocity,  geo- 
strophlc  flow  penetrated  all  the  way  to  the  bottom  of  the  equatorial 
ocean  whereas  Ekman  currents  disappeared.  Because  the  model  did  not  in¬ 
clude  zonal  dependence  and  specifically  zonal  pressure  gradients,  no 
EUC-like  features  developed.  On  the  other  hand  he  derived  an  expression 
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which  measures  the  relative  Importance  of  stratification  versus  rotation. 
For  flow  with  a  given  length  scale  L  and  depth  scale  D  in  a  thermally 
stratified  ocean  of  VSIsSIS  frequency  (gc^)1*,  rotation  (f)  will  dominate 
when  the  S  Is  small.  Clearly  in  the  equatorial  boundary  layer  where 
f-K),  stratification  takes  on  a  greater  significance  than  elsewhere.  The 
equatorial  Rossby  radius  l.e.  the  width  of  the  equatorial  boundary  layer 
can  be  defined  as  that  point  In  latitude  y*L  where  S  is  1. 
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Philander  (1973a, b)  solved  a  fully  nonlinear  set  of  momentum 
and  heat  equations  for  the  equatorial  ocean.  Though  there  were  no 
eastern  or  western  boundaries,  a  zonal  pressure  gradient  was  generated 
by  Imposing  an  eastward  sea  surface  temperature  gradient  like  that  ob¬ 
served  In  Pacific  (Colin  et  at.,  1971).  He  removed  the  zonal  dependence 
by  a  similarity  transform, then  solved  for  the  resulting  two  dimensional 
flow  by  numerical  Integration.  Results  showed  spreading  Isotherms  In  the 
vicinity  of  a  subsurface  geostrophlc  EUC  that  flowed  below  a  westward  sur¬ 
face  wind  drift  and  above  a  westward  EIC.  Inertial  terms  were  Important 
at  nearly  all  depths  on  the  equator  and  the  primary  dissipative  mechanism 
was  horizontal  friction.  McCreary's  (1980)  analysis  showed  a  similar 
vertical  current  structure.  In  addition,  flanking  the  deep  westward 
flow  and  contiguous  to  the  EUC  were  weak  tongues  of  eastward  momentum 
reminiscent  of  the  SSCC.  Unlike  Phllander's  model,  these  calculations 
were  performed  on  a  strictly  linear  set  of  equations  that  included 
vertically  diffusive  terms  as  the  only  sinks. 

Stratified  models  can  produce  geostrophic  zonal  currents  with 
realistic  shear  as  well  as  appropriate  width,  depth  and  velocity  scales 
In  the  equatorial  boundary  layer.  A  comparison  of  the  Pedlosky, 

Philander  and  McCreary  models  Indicates  that  a  baroclinic  eastward 
pressure  force  Is  necessary  to  drive  the  EUC.  By  allowing  zonal  varia¬ 
tions,  this  pressure  force  can  be  established  even  In  an  unbounded  ocean: 
McCreary  generates  It  with  a  zonally  confined  wind  stress.  Philander  by 
Imposing  a  longitudinal  surface  temperature.  McCreary's  model  would 
further  suggest  that  nonlinear  forces  can  be  neglected  as  a  first  approx¬ 
imation. 
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Through  these  efforts  we  have  discovered  some  Important  physics 
of  the  equatorial  flow  regime,  but  many  questions  remain.  We  need  to 
know  more  about  turbulent  processes  and  how  to  parameterize  them  In  large 
scale  circulation  models.  To  dramatize  this  point,  Phllander's  (1973a, b) 
similarity  solutions  predicted  an  eastward  decrease  In  width,  depth  and 
velocity  of  the  EUC.  This  may  be  a  consequence  of  his  choice  of  vertical 
eddy  coefficients  that  decrease  eastward  since  his  scaling  arguments  in¬ 
dicate  that  width,  depth  and  velocity  decrease  with  decreasing  vertical 
eddy  viscosity.  Because  Wyrtkl  and  Bennett  (1963)  argue  that  turbulent 
dissipation  Increases  west  to  east  in  the  Pacific,  It  Is  not  clear  how 
robust  the  predicted  downstream  variations  would  be  given  a  more  realis¬ 
tic  parameterization  of  diffusive  processes.  Likewise,  to  obtain  easily 
solved  equations,  McCreary  must  assume  an  eddy  viscosity  that  Is  Inversely 
proportional  to  the  VSIsSlS  frequency.  There  is  some  physical  basis  for 
this  parameterization  (Turner,  1973,  p.  143),  but  other  forms  must  be 
tested  to  gain  confidence  in  the  persistence  of  key  flow  properties. 

The  transition  from  mid-latitude  to  equatorial  dynamics  has  pre¬ 
viously  been  studied  using  matched  asymptotic  expansions  (Philander, 

1971;  Pedlosky,  1969;  Philander,  1973a, b).  This  Is  a  relatively  compli¬ 
cated  technique  that  can  sometimes  obscure  the  continuity  of  certain  flow 
features  from  one  region  to  the  next.  In  contrast,  Schneider  and  Llndzen 
(1976)  examine  thermally  driven  axl symmetric  atmospheric  circulation  over 
a  sphere  using  a  method  that  Is  uniformally  valid  at  all  latitudes.  Their 
approach  revealed  the  Importance  of  baroclinlc  pressure  gradients  in 
maintaining  well  behaved  Ekman  layers  within  the  atmospheric  equatorial 
boundary  layer.  A  uniformly  valid  representation  for  non-axlsymnetlc 
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wind-driven  flow  typical  of  the  ocean  would  certainly  be  valuable. 

In  the  following  sections  I  develop  two  analytically  tractable 
models  of  the  steady  ocean  circulation  on  a  stratified  equatorial  beta- 
plane.  Both  are  linear  and  vertically  diffusive  like  McCreary’s  model, 
but  one  dissipates  heat  by  a  Newtonian  cooling  approximation  and  the  other 
has  blharmonlc  friction.  This  work  can  also  be  viewed  as  a  generalization 
of  the  Stommel  and  Veronls  (1957)  mid-latitude  theory  to  Include  the 
equatorial  ocean.  In  that  study,  Ekman  and  geostrophlc  motions  were 
clearly  separated  Into  distinct  vertical  boundary  layers.  The  specific 
questions  I  will  address  are  1)  how  are  Ekman  and  geostrophlc  dynamics 
transformed  as  the  equator  Is  approached  and  2)  which  features  of  the 
linearized  model  persist  for  different  turbulence  parameterlzatlons.  In 
answering  these  questions  I  will  give  special  attention  to  the  role  of 
stratification  and  baroclinic  pressure  gradients. 

C.  Newtonian  Cooling  Model 

Section  C  Is  devoted  to  the  study  of  a  steady  state  model  In 
which  wind  forced  temperature  fluctuations  dissipate  via  Newtonian  cool¬ 
ing.  The  model  equations  are  developed  in  Section  Cl  (Formulation)  and 
solved  for  a  test  case  in  Section  C2  (Zonal  Winds).  Section  C3  (Variation 
of  Parameters)  explores  the  manner  in  which  results  depend  on  eddy 
viscosity,  stratification,  etc.  I  then  examine  results  for  meridional 
wind  forcing  (Section  C4)  and  conclude  with  a  brief  summary  and  discussion 
(Section  C5). 

1.  Formulation 

Consider  a  reference  atmospheric  state  of  steady,  horizontally 


uniform  downward  heat  flux  at  the  ocean  surface.  In  an  ocean  with 
horizontally  uniform  thermal  conductivity,  this  flux  will  generate  a 


background  vertical  temperature  gradient  §z  that  Is  a  function  of  depth 
only.  It  Is  hypothesized  that  deviations  from  this  atmospheric  state 

excite  perturbations  (u,v,w,p,e)  around  the  background  temperature  field 

according  to  the  linearized  set  of  equations 

-eyv  *  ^  px  -  (Au2)2 

(2a) 

syu  .  Py  .  (Av2)2 

(2b) 

P2  «  Poga0 

(2c) 

ux  +  vy  +  wz  *  0 

(2d) 

wez  =  (K0Z)Z 

(2e) 

The  coordinate  system  is  right-handed  and  subscripts  (x,y,z)  denote 
differentiation  with  respect  to  (latitude,  longitude,  depth).  I  have 
made  the  standard  equatorial  6-plane  approximation,  viz.  f=6y.  The  ocean 
Is  hydrostatic  (2c)  and  density  has  been  derived  from  an  equation  of 
state  linearized  about  temperature.  Eddy  viscosity  and  heat  diffuslvlty 
are  A  and  K  respectively,  and  will  be  assumed  constant  unless  otherwise 
stated.  Other  notation  In  (2)  is  defined  in  the  list  of  symbols  (p.  vl). 

Terms  Involving  u,v,w,  do  not  appear  In  (2).  If  the  background 
flow  Is  assumed  geostrophic  and  negligible  at  some  great  depth,  then 
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(u,v)~0.  Moreover,  from  the  resulting  vertical  advective-dlffusive 
temperature  balance,  w  -  where  Dj  is  the  depth  scale  of  the  back¬ 
ground  thermal  field.  Terms  containing  w  are  thus  small  compared  to  those 
In  (2e)  provided  that  Dj  Is  large  and/or  K  is  small.  For  the  special 
case  of  Dj  *  «  (l.e.  9Z2  =  0),  w  Is  identically  zero.  Thus  the  oceanic 
background  state  couples  to  the  perturbation  field  only  through  the  term 
we2. 

Stommel  and  Veronis  (1957)  examined  (2)  on  the  mid-latitude  f 
and  8-planes.  The  extension  of  their  work  to  the  equatorial  B-plane  is 
not  straightforward  because  for  realistic  boundary  conditions,  the  equa¬ 
tions  will  not  separate  in  y-z  (Schneider  and  Lindzen,  1976).  I  thus 
change  the  form  of  (2e)  by  making  a  Newtonian  cooling  assumption  in  this 
section 


w§2  *  -  0  (2e ' ) 

where  diffusion  of  temperature  anomalies  is  replaced  by  a  relaxation  to 
the  equilibrium  value  of  9.  The  quantity  D  is  a  depth  scale. 

I  consider  only  the  case  of  a  zonal ly  unbounded  ocean  in  which  the 
perturbation  fields  are  subject  to  the  conditions 

P0Auz  3  t^,  p0Av2  3  w  3  0  @  z  3  0  (3a) 

u  3  v  *  w  3  0  0  z  3  -  H  (3b) 


u,v,w,p,0  ■+  0  as  |y | 


-►  QO 


(3c) 
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The  depth  of  the  ocean  H  Is  taken  to  be  4000m.  Note  that  there  are  now 
no  temperature  boundary  conditions  because  temperature  follows  vertical 
velocity  ( 2e ' ) . 

To  produce  a  set  of  scaled  equations  in  which  all  terms  are  0(1), 
define  the  nondl mens Iona!  (primed)  variables 


(x,y)  *  (x\y')L,  L 


ft2gq e2 


V/« 


03K 


2  =  z'D, 


A 

7* 

A3K 

IT 

J 

82ga@_ 

z; 

(u,v,w)  «  U(u\v\£  w'),  u  =  T 


(4) 


0  ■  T  8* , 


.gaD2 


T  =  To  T* 

The  quantity  L  Is  a  barocllnlc  Rossby  radius  based  on  the  depth  scale 

A  V? 

D  =  (-gj-)  which  Is  the  Ekman  layer  depth  at  y*L*  It  can  be  derived 
from  (1)  by  setting  that  ratio  to  unity  and  A*K,  faBL.  Substituting  (4) 
into  (2)  and  dropping  primes 


•yv  +  px 


‘zz 


(5a) 


VU  +  D  *v 

jru  Ky  »2Z 


(5b) 
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ux  +  vy  +  wz  *  0  (5c) 

P2  *  e  (5d) 

w  ■  -  0  (5e) 

The  boundary  conditions  are 

uz  a  vz  a  w  =  0  @  z  ■  0  (6a) 

u  a  v  a  w  *  0  @  z  =  -  H  (6b) 

u,v,w,p,8  -*-0  as  |y|  -*•  ®  (6c) 


where  H  is  now  nondimensional  ocean  depth. 

From  (5),  a  single  expression  for  v  is  available  which  is  sixth 
order  in  the  vertical  derivative 

v(6z)  +  vx  vxxzz  ”  (vyy  "  ^  vlzz  ® 
and  which  has  separable  solutions  of  the  form 

v(x,y,z)  ■  f  S  vn(t,y,z)eUxdt 


„  $  m_  Al)z 

\  ’  jS.  Vj  U)e  ’S  *<>W 
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1  £x 

The  zonal  dependence  Is  e  because  the  ocean  Is  longitudinally  un¬ 
bounded.  Imposition  of  boundary  conditions  in  the  zonal  direction  leads 
to  a  nonseparable  partial  differential  equation  in  x-z. 

The  Hermite  functions  ipn  are  a  complete  orthonormal  set  in  lati¬ 
tude  governed  by 


(*n)yy  “  y2*n  *  "  (2n+1^n  (9a) 

.y2 

\  «  ( 2nn ! it*4) Hn(y)e  /2 

(9b) 


Hn(y)  =  (-1)"  ey2^e-y2 
n  dyn 

They  are  symmetric  (antisymmetric)  about  y=0  for  even  (odd)  n  and  decay 
(oscillate)  poleward  (equatorward)  of  a  turning  latitude  defined  from 
(9a)  as 


y  a  yt  *  ±  (2n+l)ls  (10) 

These  are  the  same  functions  that  describe  the  meridional  dependence  of 
equatorial  6-plane  waves  (Moore  and  Philander,  1977).  The  depth  depen¬ 
dence  4>*emz  obeys 


♦(6z)  +  (2n+1+*2)  <PZZ  +  ’  0 


(ID 


provided  that 


m6+  (2n+l+i2)ra2  +  U  =  0 


(12) 


Is  satisfied  for  all  n  >  0.  The  subscripts  j  in  (8)  refers  to  one  of  the 
six  roots  of  (12),  This  relation  (12)  is  similar  to  the  dispersion  of 
equatorial  waves  and  will  therefore  be  referred  to  as  a  dispersion  rela¬ 
tion.  Further  discussion  of  the  wave  analogy  follows  (21). 

From  (5),  (8)  and  the  recursion  formulae 


one  finds 


(Vy  s  ‘  Vl  +  (I)JS  Vl 


(13a) 


y*n  -  lr?-)h  Vl  +  Vl 


(13b) 


v  -  ° 


(13c) 


r 

f  ' 

/V 

un 

-l 

A 

wn 

6 

m„  j 

n,j 

A 

pn 

=  E  a  . 
J-i  n’J 

-1 

A 

0n 

-rn  . 

n,j 

>  < 

, 

Ki 


u>"  (Sr,1S 


’  * 

-1 

— m  • 
n,j 

1 

K.i  *  <!>’’  Vl 

m  . 
n,j 

. 

_ 

mn  ,z 
n,j 


(14) 


For  n  *  0  all  the  fields  (8)  and  (14)  decay  exponentially  at  high  lati¬ 
tudes  except  for  two  growing  components  of  uQ  that  are  specified  by  m2=-U. 
For  these  two  roots  m,  a  .  s  0  necessarily  and  the  dispersion  reduces 

°  I J 

to 


m4-  ii*»2  +  1=0  (15) 

There  is  a  solution  to  (5)  not  contained  in  (7),  viz.  for  v=  0 

(16a) 

(16b) 

(16c) 

(16d) 

(16e) 

Expressions  analogous  to  (14)  are 


p  =  u 
Kx  zz 


yu  +  py  =  o 


ux  +  wz  =  o 


pz  -  e 


w  =  -  e 


’  * 

’ 

A 

U  , 

1 

“1 

A 

W-1 

2 

A 

P-1 

*  ^  a-l  i 

j=i 

1 

A 

6-l 

m-l  J 

.  . 

Vy) 


(17) 
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with  dispersion 


m2  =  U 


(18) 


This  is  the  gravest  mode  in  u  and  is  designated  n=-l  since  for  that  value 
(12)  reduces  to  (18).  It  is  interesting  to  note  that  (18)  is  the  dis¬ 
persion  for  a  non-rotating  ocean  (i.e.  for  f  *  6  *  0). 

In  the  special  case  £=0,  parts  of  (8),  (14)  and  (17)  that  corres¬ 
pond  to  m=0  must  be  replaced  by 


vn_=0 


A 

un 

a„  .z  +  a„  , 
n,l  n,2 

A 

wn 

00 

-c  , 
n,l 

*  i 

Pn 

n=o 

Cn,lz  +  cn,2 

A 

en 

»  4 

Cn,l 

,  - 

where  an  and  cn  can  be  related  through  the  meridional  geostrophic  balance 
(5b). 

The  total  solution  to  (5)  is  now 

q  « Hi.!  +  E  qn]  e  UX  d£  (20) 

•'-«>  1  n*o 

A 

where  q  is  any  one  of  u,v,  etc  and  where  v_^=  0.  Provided  that 
depth  integrated  quantities 
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r-  i  ' ■  -  -> ■ 


Q 


q  dz 


(21) 


follow  readily  from  (20)  by  substituting  a„  .  -*■  a_  ./mn  4  and  subtracting 
the  series  evaluated  at  z2  from  that  at  Zj.  Depth- integrated  quantities 
will  be  discussed  in  Section  C2c. 

There  is  a  useful  analogy  between  the  present  steady  state 
frictional  model  and  time  dependent  inviscid  wave  theory.  If  w  and  e 
are  eliminated  from  (5)  and  the  time  derivative  ^  replaces  -(^j)  ,  the 
resulting  equations  are  formally  identical  to  those  that  govern  baroclinic 
equatorial  waves (e.g.  Moore  and  Philander,  1977).  This  identity  is  also 
evident  in  Figure  3  which  is  a  dispersion  plot  of  Im  (m2)  vs.  £  from  (12). 
In  wave  theory  the  ordinate  is  frequency  rather  than  vertical  wavenumber 
squared.  Though  steady  and  viscous,  the  two  families  of  curves  (n=l,2, 
3...)  designated  Ekman  and  geostrophic  are  analogous  to  inertio-gravity 
and  quasi -geostrophic  Rossby  waves  whereas  the  n=0  is  a  hybrid  equivalent 
of  the  mixed  Rossby-gravity  wave.  The  Ekman  roots  at  large  vertical 
wavenumber  will  be  strongly  surface  trapped  relative  to  the  low  vertical 
wavenumber  geostrophic  modes.  There  is  also  a  viscous  Kelvin  mode  (n=-l) 
corresponding  to  the  Kelvin  wave. 

In  addition  to  the  wave  analogy,  the  connection  with  mid-latitude 
steady  state  theories  can  be  seen  by  taking  the  limit  fc+0  in  (12): 


m"  =  -  (2n+l) 


(22a) 


m2  *  -  i£/(2n+l) 


(22b) 


Wavelength  (km) 
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Figure  3.  A  plot  of  Im(m2)  ye.  £  from  the  dispersion  relation  (12)  for  n=-l,0,l, 
2,3.  Oimensional  values  of  wavelength  (=2if£_1)  and  vertical  trapping  scale 
Re  nf 1  are  computed  for  A=K=1  cm2  sec"1,  6=1 O'"  °C  cm"1. 
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Substituting  the  expression  for  turning  latitude  (10)  Into  these  results 
In 


*■  V 

(23a) 

-  U/yt! 

(23b) 

which  is  the  dispersion  for  Ekman  (23a)  and  geostrophlc  (23b)  dynamics 
on  a  mid-latitude  6-plane  centered  at  y  a  yt  (see  Section  C2c).  An  in¬ 
teresting  feature  of  the  equatorial  vis  a  vis  mid-latitude  B-plane  Ekman 
roots  is  that  the  Ekman  depth  0g  =  Re  m"1  calculated  from  the  former  is 
not  singular  at  the  equator.  The  Kelvin  mode  is  unique  in  that  it  has  no 
mid-latitude  counterpart,  though  it  may  be  considered  an  extension  of  the 
geostrophic  roots  to  n  *  -  1  in  the  limit 

2.  Zonal  Winds 

a.  Method  of  Solution 

We  obtain  solutions  by  determining  the  unknown  amplitudes  a„  i  in 

n*j 

(8),  (14)  and  (17).  Series  expansions  for  the  boundary  conditions  are 
truncated  at  some  n=N  <»  and  values  of  u2,  vz,  w  at  za0  and  u,v,w  at 
the  ocean  bottom  are  projected  onto  them.  For  zonal  winds  characterized 
by  a  single  Fourier  component 

x  a  e  Uox  j(y)  (24) 

we  obtain  a  matrix  equation  of  the  form 


-t  ■* 

C  A  »  8 


(25) 
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whet  a  A  is  the  vector  tan  j},  8  «  (bn>  is  the  vector  projection  of  x(y) 
onto  ^n(y)  defined  by 

A  N  f  <»  ^ 

x(y)  -  Z  bni|»n(y)  ;  b  *  t(y)^  (y)dy  (26) 

n=o 

i  * 

and  C  is  a  banded  matrix  of  coefficients.  Provided  the  dimension  of  C 

is  0(100)  or  less,  (25)  can  be  solved  efficiently  using  standard  matrix 

inversion  techniques. 

We  will  examine  solutions  to  (5)  for  the  depth  invariant  quanti¬ 
ties  A=K=1  cm2  sec-1  and  i2  *  10_l*  °C  cm-1  such  that  scales  in  (4)  become 
0=1 Om,  L=46km,  U=97cm  sec-1,  W=0.21cm  sec-1,  P=4.7xl02  dyne  cm-2, 

T=2.0°C.  Zonal  wavelengths  and  vertical  trapping  scales  (Re  m-1)  then 
have  the  dimensional  values  shown  in  Figure  3.  Winds  will  be  of  the  form 

(*)  Uox 

t  =  t  *  (x,y)  =  xQe  e  (27) 


where  y  *  2xl0"3  corresponds  to  an  e-folding  distance  of  10°  latitude  and 
tq  =  0.1  dyne  cm"2.  The  value  &0  =  1.4  xlO"2  corresponds  to  a  20000km 
zonal  wavelength  which  approximates  some  of  the  largest  Trade  Wind  scales. 
A  plot  of  xv  '  is  shown  in  Figure  4.  Because  this  wind  is  symmetric  in 
latltutde,  u,w,p,6  will  also  be  even  and  v  will  be  odd.  The  vector  8 
readily  evaluated  from  (26)  is 


(28a) 


WIND  STRESS 
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Only  one  cycle  In  longitude  is  shown. 
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bn  -  0,  n  odd  (28c) 

A  small  wind  stress  curl  (y*o)  speeds  up  convergence  of  the  sequence 

$ 

therefore  minimizing  the  dimension  of  c. 

A  trial  calculation  Indicates  that  bottom  trapped  modes  (Re  m  <  0) 
are  uniformly  three  orders  of  magnitude  weaker  than  surface  trapped  modes 
for  the  winds  (27)  and  they  are  hereafter  ignored.  Series  expansions  for 
the  dependent  variables  converge  fastest  near  the  equator  and  slowest  for 
the  vertical  velocity.  The  truncation  index  is  N=77  and  allows  determi¬ 
nation  of  all  fields  to  within  0(1%).  Unless  otherwise  specified,  flow 
variables  have  been  calculated  at  the  longitude  of  maximum  westward  winds 
(x=7r  in  Figure  4). 

b.  Flow  Description 

Figure  5  shows  depth  profiles  of  various  fields  on  the  equator. 
Salient  features  include  an  eastward  EUC  (i^^  ~  120  cm  sec-1)  centered 
near  220m  below  a  surface  westward  wind  drift,  a  baroclinic  zonal  pressure 
force  which  Is  eastward  at  the  EUC  velocity  core  and  upwelling  (downwelling) 
above  (below  and  In)  this  core  correlated  with  anomalously  cold  (warm) 
water.  A  westward  EIC  below  the  EUC  appears  weak  (~5  cm  sec"1)  when  plotted 
on  the  scale  of  the  surface  current.  The  latter  is  too  strong  because 
exact  solution  does  not  permit  a  depth  variable  mixing  coefficient  which 
can  Increase  near  the  surface  in  agreement  with  the  observations  of  Jones 
(1973)  and  Crawford  (1976).  (See  however  Section  C3c).  The  pressure 
gradient  reversal  near  350m  has  never  been  observed  but  Arthur  (1960) 


Figure  5.  Depth  profiles  of  perturbation  variables  on  the  equator.  Quantities 
have  been  scaled  by  values  to  their  right  for  plotting  on  the  interval  (-1.0, 1.0). 


deduced  It  by  showing  that  the  EUC  velocity  could  not  vanish  at  depth 
otherwise. 

These  profiles  should  be  compared  to  the  flow  at  y«3.5,  poleward 
of  the  equatorial  boundary  layer.  In  Figure  6  one  sees  a  shallow  surface 
Ekman  spiral  and  a  net  meridional  transport  in  the  upper  40m.  Surface 
velocity  is  '55°  rather  than  45°  to  the  right  of  the  wind  because  of  the 
importance  of  a  geostrophic  component  of  zonal  velocity.  The  deep 
baroclinic  geostrophic  flow  is  less  strongly  surface  trapped  than  at  the 
equator  and  overall  velocities  are  weaker. 

Figure  7a  Is  a  meridional  section  of  zonal  velocity  down 
to  500m.  Only  flow  from  0  <  y  <  4  Is  shown  since  the  field  is  symmetric 
about  the  equator.  The  meridional  extent  of  the  equatorial  boundary 
layer  is  0(1  Rossby  radius)  or  about  75km  which  is  narrower  than  the  ob¬ 
served  half  width  of  the  EUC  by  a  factor  of  2.  Figure  7b  shows  the 
meridional  circulation  on  the  same  plane.  Upwelled  water  above  the  Under¬ 
current  Is  carried  poleward  in  a  surface  Ekman  layer.  Water  downwelled 
from  the  Ekman  layer  drifts  slowly  equatorward  at  depth  in  geostrophic 
response  to  the  eastward  pressure  force.  This  mass  convergence  near  the 
equator  contributes  to  both  upwelling  and  downwelling  in  the  vicinity 
of  the  EUC.  Such  a  flow  pattern  is  similar  to  that  discussed  by 
Cromwell  (1953)  and  Fofonoff  and  Montgomery  (1955),  but  because  meridional 
velocities  are  significantly  weaker  than  zonal  velocities  (especially  be¬ 
low  the  Ekman  layer)  we  might  expect  zonal  variations  to  play  an  important 
role  in  the  mass  and  momentum  balances.  Figure  8,  which  Is  a  zonal  section 
down  to  500m  of  zonal  velocity  over  one  wavelength  of  the  wind,  supports 
this  view.  At  the  surface,  u  4  0  when  ■ 


0  Implying  that  nonlocal 


Figure  7a.  Meridional  section  of  zonal  velocity.  Contour  Inter 
val  Is  40  cm  sec"1;  regions  of  westward  flow  are  stippled. 


Longitude  (rod) 


Figure  8.  Zonal  velocity  on  the  equatorial  plane  (y=0)  over 
one  wavelength  of  the  winds.  Dashed  lines  trace  velocity 
extrema  with  depth. 
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wind  effects  are  Important.  Below,  the  axis  of  maximum  velocity  (dashed 
line)  first  tilts  eastward  with  depth  down  to  220m  such  that  maximum  of 
the  eastward  EUC  lies  .2tt  radians  or  2000km  to  the  east  of  the  overlying 
westward  wind  maximum.  McCreary  (1980)  finds  a  similar  displacement  in 
his  results.  Phase  lines  then  reverse  slope  and  shift  westward  with 
depth.  I  will  elaborate  on  the  significance  of  these  results  in  the 
following  section. 

Figure  9  is  a  meridional  section  of  zonal  velocity  over  an  ex¬ 
panded  depth  range  of  200-1 500m.  The  EIC  appears  as  a  narrow  band  of 
westward  flow  below  the  swifter  EUC.  Also  evident  is  a  deep  poleward 
slanting  geostrophic  extension  of  the  EUC  suggestive  of  the  SSCC  in 
Figure  1.  Such  a  feature  appears  in  both  McCreary's  (1980)  analytical 
and  Philander  and  Pacanowski's  (1980)  numerical  calculations. 

c.  Dynamics 

Calculations  of  depth  integrated  quantities  Q  in  (21)  for  Zj=0 
arid  z2=-«  indicate  that  the  model  ocean  is  everywhere  in  Sverdrup  (1947) 
balance.  In  particular 


V  -  curl  r 

(29a) 

U  *  i£o-1(curl  t) 

(29b) 

x s  +  y(cur^ T) 

(29c) 

>y  *  -  Uo-1  y(curl  t )y 

(29d) 

This  overall  balance  masks  essential  differences  between  the  equatorial 
and  extra-equatorial  flow  field  though.  For  example,  at  mid-latitudes 
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one  can  decompose  the  Sverdrup  balance  Into  coupled  Ekman  and  geostrophlc 
components  (Stonmel,  1956).  But  as  the  equator  Is  approached,  these  com¬ 
ponents  develop  singularities  implying  that  the  internal  dynamics  must 
be  radically  transformed.  This  transformation  is  the  subject  of  the 
next  three  subsections. 

1)  Ekman  theory 

Retaining  the  equatorial  6-plane  formulation  but  requiring  that 
faBy  be  locally  constant,  one  can  rederive  Ekman' s  (1905)  theory  for 
ocean  currents  at  mid-latitudes  in  the  absence  of  horizontal  pressure 
gradients.  In  nondimens ional  notation 


-yu  =  uzz  (30a) 

yu  =  vzz  (30b) 

This  balance  yields  surface  currents  45°  to  the  right  (left)  of  the  wind 
in  the  northern  (southern)  hemisphere 


T(y)+T(x)  T(y)_x(x) 

u  = - —  ,  v  = - 7— 

(2yp  (2y  r 

and  transports  90°  to  the  winds 


Looking  for  solutions  to  (30)  of  the  form  emz,  one  obtains 


(31) 


(32) 


m4  =  -  y2  (33) 

which  Is  the  mid-latitude  equivalent  of  (23a).  The  Ekman  depth,  or  the 
depth  at  which  horizontal  velocities  are  down  to  e-1  of  the  surface 
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value.  Is  defined  from  (34)  as 

De  »  Re  m-1  ■  (p-)**  (34) 

Finally,  Ekman  pumping,  i.e.  the  vertical  velocity  at  the  base  of  the 
Ekman  layer  found  by  integrating  the  continuity  equation  Is 

*e  *  ux  +  vy  *  curl  (y)  (35a) 

where  U,  V  are  given  by  (32).  For  the  winds  (27)  this  reduces  to 

w  ■  ~2 —  +  2yt^  (35b) 

e  y 

It  is  evident  that  the  mid-latitude  theory  (30)  Is  valid  only 
for  |y|  »  0  otherwise  singularities  develop  In  (31)  -  (35).  This  Is  in 
contrast  to  the  Ekman  part  of  the  equatorial  8-plane  theory  plotted  as 
u,v  contours  on  a  meridional  plane  in  Figure  10.  Poleward  of  y»2,  solu¬ 
tions  agree  with  those  predicted  by  (30):  surface  velocity  is  45°  to 
the  right  of  the  zonal  wind,  transport  Is  all  meridional  (since  =  0) 
and  the  Ekman  depth  decreases  with  Increasing  latitude.  However,  as 
y+0,  velocities  decrease  and  become  more  meridional,  transport  decreases 
and  the  Ekman  layer  thickness  does  not  diverge. 

To  understand  how  the  equatorial  Ekman  layer  can  be  well  behaved, 

I  have  plotted  the  balance  of  forces  at  the  surface  In  Figure  11  for 
just  the  Ekman  part  of  the  total  solution.  In  the  zonal  direction,  rota¬ 
tion  balances  diffusion  at  all  latitudes  since  zonal  pressure  gradients 
associated  with  Ekman  roots  are  uniformly  small.  On  the  other  hand,  as 
y*0,  the  meridional  balance  (30b)  shifts  to 


py  ~  vzz 


(36) 


1.5 


Lotitude  (deg) 


Latitude  (deg) 


Figure  10.  Meridional  sections  of  (a)  zonal  and  (b) 
meridional  Ekroan  velocity.  Contour  Interval  Is  5  cm 
sec-1;  regions  of  westward  and  southward  flow  are  stippled. 


EKMAN  ZONAL 


EKMAN  MERIDIONAL 


(p,*vM)  (y  «■*«> 

Figure  11.  Balance  of  forces  at  the  sea  surface  computed 
from  just  the  Ekman  modes  in  the  (a)  zonal  and  (b)  meridional 
direction.  Forces  are  in  nondimensional  units. 
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In  a  constant  density  ocean,  the  balance  (36)  Is  not  possible  and  singu¬ 
larities  develop.  In  the  present  model  however,  stratification  removes 
these  by  allowing  for  barocllnlc  pressure  gradients  as  In  the  atmospheric 
model  of  Schneider  and  Llndzen  (1976).  This  result  establishes  more 
precisely  the  Importance  of  stratification  In  the  equatorial  boundary 
layer  as  predicted  by  (1). 

11)  Geos trophy 

Again  assuming  that  f  =  gy  Is  locally  constant,  one  can  calculate 
the  geostrophlc  flow  below  the  mid-latitude  Ekman  layer  from  the  non- 
dimensional  set  of  equations 


-yu  +  Px  =  0  (37a) 

yu  +  Py  -  0  (37b) 

ux  +  vy  +  wz  =  0  (37c) 

Pz  =  e  (37d) 

w  =  -  6  (37e) 


This  flow  Is  not  driven  directly  by  the  winds  but  Indirectly  via 
Ekman  pumping  (35).  Dispersion  for  e1s,°xemz  dependence  Is  given  by 

m2  .  .  uo/yz  (38) 

which  Is  the  mid-latitude  equivalent  of  (22b). 

As  with  mid-latitude  Ekman  theory,  the  balance  (37)  In  the  deep 
ocean  breaks  down  as  y*0.  For  example,  any  pressure  gradient  will  gen¬ 
erate  Infinite  velocities.  Also  the  depth  scale 
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Dg  *  fy‘\)h  (39) 

vanishes  at  y=0.  Such  singularities  do  not  appear  In  the  deep  flow  cal¬ 
culated  from  the  equatorial  8-plane  model  as  shown  In  Figure  9.  There  Is 
a  tendency  for  velocities  to  increase  and  depth  scales  to  decrease  from 
y*4  to  y*2  In  accordance  with  (37)  and  (39),  but  the  equatorial  flow  Is 
well  behaved. 

Figure  12  shows  the  balance  of  horizonal  forces  for  the  total 
solution  minus  the  Ekman  part.  For  y  >  2,  mid-latitude  geostrophy  (37) 
holds.  As  y-+0,  geostrophy  still  holds  in  the  meridional  direction 
(Figure  12b)  but  the  zonal  momentum  balance  (Figure  12a)  shifts  to 

Px  -  u22  (40) 

This  is  essentially  a  Kelvin  balance  (16).  While  the  sum  of  the  geo- 
strophic  modes  in  (20)  takes  on  Kelvin-like  characteristics  near  the 
equator  e.g.  v~0  and  u,w,p,0-4o,  by  far  the  single  most  Important  mode 
In  the  transition  to  (40)  is  the  Kelvin  mode  itself. 

Ill)  The  Kelvin  mode 

Figure  13  is  a  plot  of  profiles  at  the  equator  constructed  from 
just  the  Kelvin  mode.  Comparison  with  Figure  5  indicates  that  this  mode 
predicts  both  the  magnitude  and  vertical  structure  of  w  (away  from  the 

-y2/? 

surface),  u  and  px>  Its  Gaussian  meridional  structure  i^-e  c  accounts 
for  the  equatorial  boundary  layer  width  of  0  (1  Rossby  radius).  In 
addition  the  near  surface  eastward  displacement  of  velocity  extrema  rela¬ 
tive  to  the  winds  In  Figure  8  Is  due  to  the  Kelvin  phase  relationship 

cxz“  a!  *  ■ Im  hr*  *  *  (TJ  (41) 


TOTAL  minus  EKMAN  ••  MERIDIONAL 


Latitude 


Latitude 


Figure  12.  Balance  of  forces  at  the  sea  surface  computed  from 
all  modes  excluding  the  Ekman  modes  In  the  (a)  meridional  and 
(b)  zonal  direction.  Forces  are  In  nondlmenslonal  units. 


Figure  13.  Depth  profiles  of  perturbation  variables  on  the  equator  computed  from 
just  the  Kelvin  mode. 


which  Is  found  by  setting  Im(mz  +  UQx)  *  constant  and  using  (18).  (The 
westward  tilt  below  250m  Is  due  to  geostrophic  components  that  according 
to  (22b)  extend  deeper  and  have  Cx2>0) .  In  contrast  to  constant  density 
models  which  ah  initio  set  ux*0,  (41)  reflects  the  fact  that  ux  Is  Im¬ 
portant  and  that  =  0  In  the  Kelvin  continuity  balance  (16c).  The 
physical  Implication  of  such  zonal  variability  is  that  flow  at  a  particu¬ 
lar  depth  depends  not  only  on  local  surface  conditions  but  also  conditions 
at  upstream  and  downstream  positions. 

At  mid- latitudes,  winds  project  entirely  onto  the  surface  Ekman 
layer  which  in  turn  drives  the  deep  ocean  via  Ekman  pumping  as  described 
above.  In  contrast,  zonal  wind  stress  forces  the  deep  equatorial  ocean 
directly  through  projection  onto  Kelvin  and  geostrophic  modes  (Figure  14). 
The  equatorial  Ekman  layer  is  thus  of  secondary  importance  in  transferring 
zonal  momentum  from  the  wind  field  to  the  ocean  which  accounts  for  the 
weak  horizontal  currents  near  the  equator  in  Figure  10.  Its  primary  role 
Is  to  balance  a  vertical  velocity  at  the  surface  associated  with  the 
Kelvin  (Figure  13)  and  weaker  geostrophic  modes. 

The  Kelvin  balance  (16)  is  not  possible  in  a  constant  density 
ocean.  Key  features  of  the  equatorial  circulation  related  to  this  mode 
e.g.  geostrophic  balance  of  vertically  sheared  zonal  currents,  width  of 
the  boundary  layer,  etc.  depend  on  stratification.  These  results  confirm 
the  expectation  from  (1)  that  scratifi cation  should  be  important  not 
only  in  the  equatorial  Ekman  layer  but  also  In  the  deeper  equatorial 
ocean  as  wel 1 . 

3.  Variation  of  Parameters 

So  far  we  have  considered  solutions  only  for  a  single  zonal 
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wavenumber  l0  corresponding  to  a  wavelength  of  2x10“  km  In  the  winds. 
Since  observed  winds  consist  of  a  spectrum  of  wavenumbers,  it  is  of 
interest  to  know  how  sensitive  the  results  are  to  l0  and  what  contribu¬ 
tion  other  wind  scales  will  make.  Furthermore,  separation  of  variables 
in  the  model  equations  requires  that  mean  stratification  and  diffusivities 
be  uniform  with  depth.  In  the  ocean  however,  these  quantities  can  vary 
over  several  orders  of  magnitude.  A  typical  temperature  profile  would 
show  a  surface  mixed  layer  (6Z~0)  tens  of  meters  thick  above  a  sharp 
thermocline  (§Z~10"3°C  cm-1);  deeper  water  is  less  stratified 
(5Z~10“  -10“*  °C  cm-1).  Gregg  (1976)  using  a  Cox  number  model,  calcu¬ 
lated  eddy  heat  diffusivities  between  0.15-5.0  cm2  sec-1  at  155°W  on  the 
equator.  Likewise  Robinson  (1966),  Jones  (1969)  and  Crawford  (1976) 
using  empirical  formulae,  found  A  to  range  between  -1-100  cm2  sec-1.  In 
addition,  these  studies  indicate  that  the  turbulent  Prendtl  number  Pr=K/A 
is  less  than  unity  consistent  with  the  notion  that  stratification  inhibits 
the  transfer  of  heat  more  effectively  than  momentum  (Turner,  1973). 

In  this  section  I  investigate  how  depth,  width,  velocity  and 
other  properties  of  the  equatorial  circulation  vary  with  these  parameters. 
The  control  case  is  that  of  zonal  wind  forcing  discussed  in  Section  2; 
the  procedure  will  be  to  change  only  one  of  A,Pr,ez  or  lQ  holding  all 
the  rest  fixed.  The  section  concludes  with  a  boundary  layer  calculation 
in  which  dlffusivity  and  stratification  vary  continuously  with  depth  on 
the  equatorial  plane. 

a.  Zonal  Wavenumber  (tQ) 

Because  near  the  equator  zonal  winds  project  directly  onto  the 
Kelvin  mode,  one  might  expect  from  the  dispersion  relation  (18)  that 
increasing  the  zonal  wavelength  increases  the  penetration  of  surface 


effects.  One  can  see  these  developments  In  Figure  15  in  which  is  a  plot 
of  zonal  velocity  ye.  depth  at  the  equator  for  various  lQ.  Solid  curves 
were  calculated  by  the  method  outlined  in  Section  C2a;  the  dashed  curve, 
u  =  (z  +  H),  is  a  direct  calculation  of  the  zonal  momentum  (Auzz=0) 
for  the  special  case  iQ  =  0.  Numbers  in  parentheses  are  wavelength  in 
thousands  of  kilometers. 

Over  two  orders  of  magnitude,  lQ  ~  2lT/2*103km  -  2Tr/2-<105km,  the 

depth  of  maximum  eastward  current  increases  as  £  in  accordance  with 

o 

(18).  This  trend  continues  until  currents  feel  the  bottom.  Then  further 
decrease  to  &  =  0  leads  to  a  linear  shear  profile  over  the  full  depth 
range  of  the  ocean.  Note  however  that  the  approach  to  this  limit  is 
slow:  currents  driven  by  winds  with  a  zonal  wavelength  of  2xl05km  pene¬ 
trate  to  only  about  ~2km.  Figure  15a  also  indicates  that  for  currents 
that  do  not  extend  the  bottom,  the  velocity  decreases  as  £q  .  The 
width  of  the  equatorial  boundary  layer  does  not  change  with  however 
(Figure  15b)  nor  does  the  geostrophic  balance  of  zonal  currents  below  the 
shallow  equatorial  Ekman  layer. 

In  the  present  model  in  which  there  are  no  meridional  boundaries, 
zonal  pressure  gradients  along  the  equator  are  generated  by  zonal  wind 
convergences.  Thus  changing  £q  in  the  wind  field  changes  px  *  U0p  in 
the  ocean.  One  can  then  infer  from  Figure  15  that  baroclinic  zonal 
pressure  gradients  are  important  a)  in  maintaining  subsurface  structure 
like  the  EUC,  EIC,  SSCC,  etc.  and  b)  in  limiting  the  strength  of  zonal 
flow.  McCreary  (1980)  arrived  at  similar  conclusions. 

b.  Eddy  Dlffusivities  (A,K)  and  Stratification  (9Z) 

Table  1  summarizes  the  dependence  of  various  scales  on  eddy 
dlffusivities  and  stratification.  One  point  to  emerge  is  that  increasing 


ZONAL  WAVENUMBER 

U (21/210  cm  sec**  — - 

U  (20)/670  "  - 

U(200)/2100  "  - 

u(®)/40,000  "  - 


-(.0  -0.5  0  0.5 


Figure  15.  Effects  of  varying  zonal  wavenumber  (£q)  on  (a)  zonal 
velocity  at  the  equator  and  (b)  width  of  the  equatorial  current 
system  plotted  as  ipQ.  Numbers  in  parentheses  are  zonal  wavelength 
(s2nZo~2)  in  thousands  of  kilometers.  The  three  curves  in  (b) 
overlap. 
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stratification  is  Indistinguishable  from  decreasing  heat  dlffuslvlty  and 
vice  versa  since  ©z  and  K  always  appear  in  the  scales  (4)  as  a  ratio 
(K/0Z)W.  A  second  point  is  that  in  the  limit  of  constant  density  0Z  -►  0, 
several  scales  become  singular.  However  when  the  stratification  is  so 
weak  i.e. 

<  («) 


such  that  the  depth  scale  D>H,  the  heat  balance  is  inconsequential.  The 
ocean  Is  at  this  stage  essentially  constant  density  and  the  appropriate 
scales  become 


D  =  H  , 


U  * 


H  T 


o’ 


ill 


PH4 

-p- 


P  = 


_  A 


(43) 


The  dependence  of  these  scales  on  A  is  shown  parenthetically  in  Table  1. 
Clearly  the  stratified  model  is  less  sensitive  to  this  parameter. 

Figure  16  is  analogous  to  Figure  15  except  that  here  we  have 
varied  A,  keeping  lQ,  0Z  and  Pr  fixed.  Numbers  in  parentheses  are 
eddy  momentum  (and  heat)  dlffusivitles  In  cm2  sec"1.  The  velocities 
attenuate  and  extend  to  greater  depths  as  turbulent  diffusion  increases 
and  the  current  system  likewise  expands  latitudinally.  As  expected  from 
Table  1  though,  the  solutions  vary  much  more  slowly  than  A.  In  particu¬ 
lar  the  width  of  the  equatorial  flow  field  Increases  by  a  factor  of  only 
2.5  when  A  Is  changed  by  two  orders  of  magnitude. 

Figure  17  shows  the  effect  of  varying  stratification  for  fixed 
A,  Pr  and  lQ.  Numbers  in  parentheses  are  vertical  temperature  gradients 
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DIFFUSI VITY 

U(1)/670  cm  sec'1  - 

U(10)/130  "  - 

U  (100)/ 27  "  - 


*1.0  -0.5  0.  0.5 


Figure  16.  Effects  of  varying  dlffuslvlty  (A,K)  on  (a)  zonal 
velocity  at  the  equator  and  (b)  width  of  the  equatorial  current 
system  plotted  as  Numbers  In  parentheses  are  dl.ffuslvltles 
In  cm2  sec-1. 
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Figure  17.  Effects  of  varying  stratification  (02)  or  Prandtl  number 
(Pr)  on  (a)  zonal  velocity  at  the  equator  and  (b)  width  of  the 
equatorial  current  system  plotted  as  ip  .  Numbers  In  parentheses 
are  either  vertical  temperature  gradient  In  °C  per  100m  or  Pr*1. 


In  °C  per  100  m.  The  dashed  line,  u  *  Sp’TF  +  H)2»  ^as  ^een  ca^cu^ate<l 
directly  from  the  zonal  momentum  equation  for  the  limiting  case  0Z  ■  0 
on  the  assumption  that  «  t^(0).  One  can  see  that  Increasing 

stratification  limits  both  the  depth  and  strength  of  the  zonal  flow,  but 
the  variation  with  §z  Is  weak  as  In  the  case  of  eddy  momentum  dlffuslvlty. 
In  fact  to  approach  the  limit  case,  the  basic  stratification  must  satisfy 
8Z  <  10"n  °C  cm-1  for  A=Kalcm2  sec-1  and  H=4000m.  As  seen  In  Figure  15b 
stronger  stratification  also  leads  to  an  Increase  In  the  width  of  the 
flow,  but  again  the  dependence  is  very  weak.  The  width  corresponlng  to 
@z  =  0  Is  not  plotted,  according  to  (43)  it  Is  only  about  25  cm!  Further¬ 
more,  in  this  limit  zonal  currents  at  the  equator  are  not  geostrophic 
as  discussed  in  the  previous  section. 

As  mentioned  earlier,  the  effects  of  changing  stratification  are 
Indistinguishable  from  changing  the  heat  diffusivity.  Thus  one  can 
alternately  interpret  Figure  17  in  terms  of  varying  the  Prandtl  number 
for  fixed  A,  ©z  and  «,Q.  The  numbers  of  parentheses  are  then  simply 
Pr-i. 

It  is  interesting  to  note  that  while  D,  L  and  U  all  vary  with  A, 

K  and  0Z,  the  scale  for  mass  transport  UDL  (*t0/p08)  Is  independent  of 
these  parameters. 

c.  An  Approximate  Analysis 

I  now  develop  a  simple  model  that  allows  both  stratification 
and  eddy  diffuslvlties  to  vary  continuously  with  depth.  In  this  way  one 
can  Introduce  a  surface  mixed  layer  into  the  problem.  The  analysis  is 
based  on  the  dominance  of  the  Kelvin  mode  near  the  equator  and  is  valid 
only  on  the  equatorial  plane  (l.e.  y*0). 
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Assuming  v  =  0, 


Px  *  (A'uz)z 

(44a) 

yu  +  py  *  0 

(44b) 

ux  +  w2  =  0 

(44c) 

Pz  -  0 

(44d) 

k* 

"■'•Ti8 

(44e) 

where  A1,  K* ,  6^  are  nondlmenslonal ,  depth  variable  quantities  scaled  by 

their  depth  averaged  dimensional  values  A,  K, 

Ad1m  *  A  A> 

(45a) 

Kd1m  *  K  K’ 

(45b) 

<®z>d1m  a 

(45c) 

If  the  vertical  velocity  boundary  conditions  are  relaxed,  (44)  can  be 

solved  at  the  equator  where  winds  project  entirely  onto  Kelvin- like  flow 

(Figure  14)  subject  to 

Auz  ■  @2  =  0 

(46a) 

u  ■  0  @  z  3  -  H 

(46b) 

Combining  (44c-e)  leads  to 

ux  "  ^z 

(47) 

which  by  comparison  with  (44a)  yields  a  well  posed  problem  If 


u  *  p  ;  A' 


K' 

si 


(48) 


For  the  winds  In  Figure  4,  (44a),  (47)  and  (48)  lead  to 

A1  l 

uzz  +  7T  uz  _  1  T  u  *  0 


(49) 


This  equation  Is  trivially  solved  by  standard  Runge-Kutta  techniques  for 
arbitrary  A'(z). 

I  analyze  two  cases  characterized  by  depth  averaged  parameter 
values  of  A*K*1cm2  sec-1  and  9Z  =  10-4  °C  cnr1  such  that  the  same  scales 
(4)  apply  to  both.  First,  assume  that  A1,  K',  @z  are  depth  Independent 
which  Is  simply  the  problem  considered  more  exactly  In  Section  C2.  The 
solution  to  (49)  Is  plotted  as  a  light  solid  line  In  Figure  18.  One  can 
see  by  comparison  with  Figure  5  that  this  curve  approximates  both  the 
amplitude  (to  within  1555)  and  vertical  structure  of  zonal  velocity.  Now 
compare  this  control  case  to  one  in  which 


A'  *  max  (0.51 

,  lOOe  6z) 

(50a) 

K'  *  1 

(50b) 

6^  =  min  (1.96 

,  0.01e'6z) 

(50c) 

where  min  (max)  Implies  the  minimum  (maximum)  value  in  parentheses  is 
taken  and  6  *  0.5.  In  dimensional  terms,  this  recipe  (50)  calls  for 
strong  mixing  near  the  surface  (A^  ~  100  cm2  sec-1)  In  a  shallow  mixed 
layer  (~10-s  °C  cm"i)  below  which  viscosity  and  stratification  are  con¬ 
stant  at  0.51  cm2  sec-1  and  1.96xl0-4  °C  cm-1.  Plots  of  A^  and  6^ 


appear  In  Figure  18. 

The  heavy  line  In  Figure  18  Is  the  solution  to  (49)  for  (50). 
Compared  to  the  control  case,  surface  velocities  are  much  reduced  and 
more  slab-like.  McCreary  (1980)  finds  similar  behavior  in  the  mixed 
layer  of  his  model.  Conversely  at  the  level  of  the  EUC  where  parameter 
values  are  comparable,  flow  is  much  less  modified.  Slower  decay  of 
current  speed  below  the  EUC  in  the  first  case  can  be  ascribed  to  the 
slightly  higher  diffusivity  there  (1cm2  sec-1  vs.  0.51  cm2  sec-1)-  In 
both  cases  there  Is  a  nonzero  vertical  velocity  at  the  surface  (not 
shown)  that  would  be  nullified  if  a  shallow  Ekman  layer  were  Included. 

4.  Meridional  Winds 

In  this  section  I  examine  the  solution  to  (5)  for  meridional 
winds  of  the  form 

T<y)(x.y)  -  T0  e  'V  e-W!  (51) 

which  is  identical  to  (27)  except  for  the  direction.  Other  physical 
parameters  are  the  same  as  in  the  zonal  wind  case.  The  method  of  solu¬ 
tion  is  that  discussed  In  Section  C2a  but  now  vz  (z=0)  is  projected  onto 
even  t/>n.  Thus  v  has  even  symmetry  in  latitude  and  u.w.p.e  have  odd 
symmetry. 

a.  Flow  Description 

Because  only  v  Is  nonzero  at  the  equator,  Figure  19  displays 
profiles  of  dynamical  variables  In  the  equatorial  boundary  layer  at 
y»0.5.  Winds  from  the  north  drive  a  surface  flow  25°  to  the  right  of 
the  forcing  (or  by  symmetry  to  the  left  In  the  southern  hemisphere)  and 
simultaneously  upwelllng  (downwelllng  in  the  south).  Below  the  south- 


ward  surface  drift  is  a  northward  return  flow.  Zonal  pressure  gradients 
are  not  important  and  there  is  no  significant  structure  below  40m.  For 
nondlmensional  latitudes  y  >  1,  flow  is  similarly  confined  to  the  upper 
40m  in  a  classical  Ekman  layer.  Figure  20  shows  the  circulation  on  a 
meridional  plane  at  the  longitude  of  maximum  southward  winds.  North 
(south)  of  the  equator,  upwelling  (downwelling)  is  forced  by  meridional 
divergence  (convergence)  whereas  at  the  equator  there  is  no  vertical 
motion.  Cromwell  (1953)  argued  for  a  similar  circulation  pattern  under 
the  influence  of  meridional  winds  based  on  various  tracer  distributions 


in  the  equatorial  Pacific. 

Aside  from  obvious  symmetry  characteristics,  two  features  dis¬ 
tinguish  the  response  of  an  equatorial  ocean  to  meridional  vs.  zonal 
winds.  First,  for  the  same  magnitude  wind  stress  t  ,  zonal  forcing 
generates  currents  an  order  of  magnitude  stronger  than  meridional  as 
seen  from  the  relative  scales  in  Figure  5  and  Figure  19.  Second,  zonal 
wind  effects  penetrate  deeper.  One  can  clarify  these  differences  by 
comparing  the  dynamics  of  each  case, 
b.  Dynamics 

As  in  the  case  of  zonal  winds,  depth  integrated  quantities  for 
forcing  obey  the  Sverdrup  relations  (29).  Moreover,  the  Ekman  com¬ 
ponent  of  flow  undergoes  a  similar  transformation  as  y-K)  viz.  -yu  =  uz2 
everywhere  but  (yu  a  vzz)  -*•  (py  =  vzz)  in  the  equatorial  boundary  layer. 
However  the  winds  (51)  do  not  excite  a  significant  geostrophic  component. 
At  mid-latitudes  this  is  due  to  weak  Ekman  pumping 


which  compared  to  wg  for  zonal  winds  (35b)  is  0(yJto)«l.  Near  the 


Figure  20.  Meridional  circulation  in  the  equatorial  boundary  layer  for  southward 
winds.  Scales  for  v  and  w  shown  at  lower  right  are  50  cm  sec"1  and  2.0x10” 3  cm 
sec”1  respectively 
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equator  zonal  winds  project  directly  onto  Kelvin-like  flow  which  is  in 
cross-stream  geostrophic  balance  and  which  extends  to  several  hundred 
meters  depth.  Conversely,  meridional  winds  project  almost  entirely  onto 
the  gravest  mode  in  v,  i.e.  n=0  in  Figure  3.  As  £o-K),  the  dispersion  (15) 
for  this  mode  does  not  yield  a  geostrophic  root  so  that  the  dynamics  are 
predominantly  frictional  and  strongly  surface  trapped.  The  frictional 
nature  of  the  response  also  accounts  for  its  weakness:  even  modest  de¬ 
formations  of  the  thermal  field  near  the  equator  if  geostrophically 
balanced  would  lead  to  much  larger  velocities. 

A  more  subtle  distinction  between  zonal  and  meridional  wind  forc¬ 
ing  is  that  in  the  latter  case,  flow  at  any  given  longitude  is  influenced 
by  local  winds  only  and  not  by  winds  at  neighboring  longitudes.  Plots  of 
variables  on  the  x-z  plane  exhibit  no  change  of  phase  with  depth  i.e. 

Cxz  =  ®  Fl9ure  8)*  This  is  due  to  the  fact  that  for  the  Ekman  roots, 
not  only  is  p  negligible  in  (5a)  but  also  u  is  negligible  in  the  con- 
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tinuity  equation  (5d). 

5.  Discussion 

The  essential  ingredients  of  the  present  steady  state  theory  are 
variable  rotation  and  continuous  density  stratification.  The  linearized 
mathematical  formulation  precludes  meridional  boundaries  and  allows  only 
for  highly  idealized  types  of  winds.  In  spite  of  its  simplicity  though, 
results  show  much  in  conroon  with  oceanic  observations.  For  a  range  of 
parameters,  realistic  width,  depth  and  velocity  scales  are  predicted 
that  do  not  depend  too  strongly  on  the  eddy  coefficients.  Geostrophic 
zonal  flows  reminiscent  of  the  SEC,  EUC,  EIC  and  SSCC  can  be  generated 
along  with  convincing  meridional  circulation  patterns.  Baroclinic 


zonal  pressure  gradients  provide  the  source  of  momentum  for  deep  zonal 
flows  near  the  equator. 

In  Section  C2c  it  was  shown  that  meridional  baroclinic  pressure 
gradients  maintain  a  well  behaved  equatorial  Ekman  layer  contrary  to 
the  qualitative  notion  stemming  from  constant  density  models  that 
singularities  should  develop  there.  In  the  case  of  zonal  winds,  the 
Ekman  layer  is  characterized  by  weak  horizontal  currents  because  most 
or  all  of  the  wind  stress  at  a  particular  latitude  projects  directly  onto 
Kelvin  and  Kelvin-like  geostrophlc  flow.  This  geostrophic  flow  diffuses 
vertically  (unlike  at  mid-latitudes)  and  is  associated  with  a  surface 
vertical  velocity  which  the  Ekman  layer  nullifies.  The  response  to 
large  scale  meridional  winds  on  the  other  hand  is  entirely  confined  to 
a  shallow  surface  frictional  layer. 

Baroclinic  meridional  pressure  gradients  that  develop  in  the 
equatorial  Ekman  layer  imply  a  temperature  distribution  that  is  not  in 
geostrophic  balance.  Thus,  shallow  temperature  anomalies  due  to 
frictional  effects  represent  a  source  of  error  in  the  geostrophic  calcu¬ 
lation  of  near  surface  ^uatorial  flows.  This  error  can  be  severe  in 
the  presence  of  winds  with  a  strong  southerly  or  northerly  component 
since  frictionally  induced  meridional  temperature  gradients  are  nonzero 
even  at  the  equator  in  this  case.  For  strictly  zonal  winds,  the  rela¬ 
tive  size  of  Ekman  and  geostrophic  pressure  gradients  (py)  shown  in 
Figures  11  and  12  indicates  an  error  of  n,10%  at  y=l. 

The  model  suggests  deep  baroclinlcity  like  that  observed  by 
Luyten  and  Swallow  (1976)  and  Luyten  and  Eriksen  (1980).  From  (22b)  and 
(38),  the  vertical  wavenumber  at  high  latitudes  Is  very  small  Implying 
•  strongly  barotropic  interior  flow.  As  the  equator  is  approached,  the 
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wavenumber  progressively  Increases  and  asymptotes  to  the  Kelvin  wavenumber 
given  by  (18)  for  which  the  vertical  wavelength  Is  only  a  few  hundred 
meters.  The  model  differs  from  the  observations  though  In  that  (18) 
predicts  a  much  more  rapid  decay  from  the  surface. 

Meridional  winds  or  zonal  winds  with  a  meridional  shear  produce 
zonal  currents  antisymmetric  around  the  equator.  In  a  constant  density 
ocean,  either  of  these  in  combination  with  easterlies  leads  to  an  EUC 
structure  centered  off  the  equator.  In  a  stratified  ocean,  the  response 
to  antisymmetric  zonal  winds  (which  has  been  calculated  but  not  presented) 
and  meridional  winds  (Section  C4)  Is  confined  to  a  shallow  surface  Ekman 
layer.  Conversely,  geostrophlc  flow  driven  by  zonal  wind  stress  extends 
much  deeper  (Section  C2).  Thus  In  the  present  model,  meridional  or  anti¬ 
symmetric  zonal  winds  do  not  result  In  an  asymmetric  EUC  structure.  Such 
structure  In  a  linear  calculation  is  then  only  possible  If  time  dependence 
Is  Introduced  (q.v.  Chapter  III).  On  the  other  hand, Phi lander  (1973b) 
suggested  that  In  a  stratified  ocean  nonlinear  effects  may  produce  a 
steady  displacement  of  the  EUC  off  the  equator  In  response  to  meridional 
winds. 

Observations  suggest  that  the  maximum  speed  of  the  EUC  lies  to 
the  east  of  the  maximum  westward  wind  stress  as  predicted  by  the  model 
(Figure  8).  Current  measurements  between  140°W  and  90°W  (Knauss,  1960) 
Indicate  a  broad  EUC  velocity  maximum  centered  near  115°W.  Long  term 
averages  of  the  Pacific  Trade  Wind  field  along  the  equator  (Knauss,  1963; 
Meyers,  1979)  show  maximum  easterlies  located  near  140°W.  This  agreement 
with  the  model  may  be  fortuitous  however  since  the  current  measurements 
were  made  only  In  Spring,  1958.  If  wave  processes  have  contaminated  the 


record,  the  pattern  of  zonal  velocity  so  described  may  not  be  representa¬ 
tive  of  the  long  term  average.  No  suitable  combination  of  velocity  and 
wind  data  are  available  In  the  Atlantic  to  compare  with  the  theory. 

Even  though  observations  Indicate  that  nonlinearity  Is  Important 
In  the  EUC,  this  and  McCreary's  (1980)  strictly  linear  calculation  pro¬ 
duce  remarkably  realistic  flow  patterns.  Hence  as  with  constant  density 
models.  Inertial  terms  derived  from  a  linear  theory  may  predict  the  effects 
of  nonlinearity.  For  the  case  of  zonal  winds  (27),  I  computed  uu  ,  vu  , 

A  Jr 

wu2  (zonal  momentum),  uvx>  Wy,  wvz  (meridional  momentum)  and  u0x,  vey, 
w©  (heat  balance)  as  a  function  of  latitude  and  depth  at  the  longitude 
of  maximum  easterlies.  Results  show  that  quadratic  terms  In  the  heat 
and  meridional  momentum  balances  are  small  relative  to  the  linear  terms. 

In  the  zonal  direction,  plots  of  nonlinear  terms  relative  to  the  dominant 
linear  term  (uzz  near  the  surface,  px  below)  reveal  that  nonlinearity  is 
Important  In  the  equatorial  boundary  layer  but  not  at  mid-latitudes. 

For  example  Figure  21  shows  that  at  the  core  of  the  EUC  below  the  west¬ 
ward  wind  maximum  there  is  both  equatorward  advection  of  eastward  momen¬ 
tum  (vuy>o)  which  tends  to  narrow  the  flow  and  downstream  advection 
(uux>o)  which  tends  to  reduce  its  strength.  Above  and  below  the  core, 
there  Is  a  gain  In  eastward  momentum  (wuz<  o)  which  tends  to  Increase 
the  vertical  extent  of  the  EUC.  A  comparison  of  linear  and  non¬ 
linear  numerical  calculations  In  Philander  and  Pacanowski  (1980, 

Figure  12)  reveals  similar  tendencies. 

The  model  shows  spreading  Isotherms  near  the  equator,  but  It 
also  develops  temperature  Inversions  (Figure  22).  The  near  surface 
ones  are  due  to  the  restriction  that  e  «  w  °  0  at  z  *>  0.  Down- 
welling  inmedlately  below  the  surface  can  then  carry  warmer  water  to 


NONLINEAR  TENDENCIES 
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Figure  21.  Nonlinear  tendencies  as  a  function  of  non- 
dimensional  latitude  at  the  depth  of  the  EUC  core. 
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Figure  22.  Meridional  section  of  total  temperature.  Contour 
Interval  Is  0.33  °C;  the  stippled  area  represents  water  warmer 
than  the  surface  temperature.  Superimposed  are  Isotachs  of 
eastward  velocity  In  the  EUC  (dashed  lines). 
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depth  raising  ambient  temperature  above  the  surface  value  which  cannot 
change.  Above  the  EUC  the  Inversions  are  due  to  stronger  vertical 
advection  of  the  background  temperature  field  near  the  surface  than 
below..  Colder  water  Is  thus  carried  from  greater  depth  to  rest  over 
warmer  water.  McCreary's  (1980)  linearized  heat  balance  leads  to  sim¬ 
ilar  inversions.  These  static  instabilities  imply  strong  mixing 
through  a  convective  adjustment  process  and  suggest  that  the  parameter¬ 
ization  of  turbulent  conduction  in  this  study  is  inadequate.  The  next 
section  explores  an  alternate  parameterization  scheme. 
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D.  Biharmonic  Friction  Model 

Under  this  heading  I  develop  a  model  which  Incorporates  the  con- 
ventlonal  Laplaclan  turbulent  heat  diffusion  (proportional  to  a  second 
derivative)  as  well  as  a  less  comnon  biharmonic  vertical  friction  (pro¬ 
portional  to  a  fourth  derivative).  Though  solutions  will  be  compared 
with  the  Newtonian  cooling  case,  the  purpose  Is  not  to  demonstrate  the 
superiority  of  one  type  of  parameterization  over  another.  It  Is  rather 
to  show  that  Important  dynamical  balances  are  Insensitive  to  details  of 
these  schemes.  The  particular  formulation  discussed  In  Section  D1  leads 
to  separable  solutions  on  the  y-z  plane  as  In  the  Newtonian  cooling  case. 
Section  D2  Is  an  examination  of  the  response  to  zonal  winds  like  those 
of  (27)  and  Section  03  discusses  thermally  forced  solutions.  Section  D4 
Is  a  brief  recapitulation  and  commentary. 

1.  Formulation 

Using  the  same  notation  as  in  Section  6,  governing  equations  are 


-eyv  +  Px  -  -  (AD2uzzz)  (53a) 

Byu  +  J_  P  =  -  (AD2VZZZ)Z  (53b) 

Pz  *  Po9a0  (53c) 

Ux  +  Vy  +  wz  0  (53d) 

wez  -  (Kez)z  (53e) 

subject  to 


-‘>0AUZZ2  ’  T<X)  •  V'zZZ  '  ■'<y>  •  6Z  -  ° 


@  z«0 


(54a) 
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u  ■  v  «  w  ■  02  ■  u22  ■  v22  ■  0  @  z»-H  (54b) 

u,v,w,p,0  0  as  |y|  -*> «  (54c) 

Blharmonlc  friction  distinguishes  this  model  from  nearly  all  other  ocean 
general  circulation  models  which  typically  use  Laplaclan  vertical  momentum 
diffusion.  The  blharmonlc  formulation  Is  less  unusual  In  parameterizing 
horizontal  turbulence  In  numerical  eddy  resolving  models  (e.g.  Holland, 
1978)  where  It  strongly  damps  grid  scale  fluctuations  but  not  mesoscale 
eddies.  The  logic  for  Its  Introduction  Into  (53)  Is  that  the  equations 
now  separate  In  all  three  directions. 

The  form  of  the  stress  boundary  conditions  In  (54a)  guarantees 
that  the  ocean  will  be  In  Sverdrup  balance.  However,  the  additional 
derivatives  In  (53a, b)  require  another  boundary  condition  on  u  and  v  at 
the  top  and  bottom.  I  have  assumed  the  null  hypothesis  for  velocity 
curvature, l.e.  u22  *  v22  *  0.  Temporarily  the  perturbation  heat  flux 
(proportional  to  ez)  is  set  to  zero  at  both  upper  and  lower  boundaries, 
though  the  upper  condition  can  be  relaxed  or  replaced  by  a  surface  tem¬ 
perature  condition  (Section  03).  Hereafter  A  and  K  will  be  assumed 
constant. 

From  (53)  and  (54)  one  can  produce  a  set  of  equations  In  which 
all  terms  are  0(1)  by  Invoking  the  scaling  arguments  (4).  A  single 
12th  order  equation  In  v  follows 

v(12z)  *  vx  -»«(4z)  -  t’yy  *  >*wJ<4z)  *  0  <55) 

which  has  separable  solutions  of  the  form 
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v(x,y,z)  *  I  S  (£,y,z)eUx  d£ 
J.oo  n*o  n 


v  mn  i(*)2 

•  -  Z  a_  <(£  e  n,J  i|»  (y) 
n  j*i  n,jv  rnv,,/ 


^  12 
V 


(56) 


The  Hermite  functions are  again  governed  by  (9)  but  the  vertical  de¬ 
pendence  0*emz  obeys 


*(12z)  +  <2n+1+JlZ)  4»(4z)  *  1W  -  0  (57) 


provided  that 


m12  +  (2n+l+£2)m4  -  i£  =  0  (58) 

holds  for  n  >  0.  The  quantities  u,w,p,0  follow  from  (56)  and  the  non- 
dimensional  version  of  (53). 

As  with  the  Newtonian  cooling  model  (Section  C),  there  are  roots 

of  (58)  for  n  =  0.1. e.  m4  «  U  that  results  in  ill  behaved  u  at  high 

o 

latitude.  Eliminating  these  leads  to  an  analogue  of  (15) 

ro*  +  l£m4  +  1  *0  (59) 

In  addition  there  Is  an  n*-l  Kelvin  mode  with  v  =  0  and 

m4  *  -  l£  (60) 

We  could  have  obtained  (55)  and  (57)- (60)  from  their  equivalents 
In  Section  Cl  by  the  transformation  +  1(^y)2,  m  -*■  im2.  Thus  the 
analogy  with  equatorial  waves  as  well  as  with  mid-latitude  Ekman  and 
geostrophic  theory  is  preserved.  Specifically  as  £  -*  0,  (58)  reduces  to 
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-  (2n+l )  =  -  yt2 

(61a) 

u/(2n+l)  -  u/yt2 

(61b) 

The  roots  (61a)  Imply  a  different  vertical  structure  than  traditional 
Ekman  theory  (23a)  but  they  are  still  surface  trapped,  steady  state, 
frictional  analogues  of  inertia-gravity  waves.  Likewise  (61b)  corresponds 
to  the  steady  state  Rossby  wave  analogues  (22a)  and  (22b). 

2.  Zonal  Winds 
a.  Method  of  Solution 

Equations  (53)  and  (54)  are  solved  according  to  the  method  out¬ 
lined  in  Section  C2a  for  winds  of  the  form  (27)  shown  in  Figure  4  and  for 
parameter  values  A=K=40cm2  sec-1,  6z“10~*  °C  cm"1.  The  values  of  A,K 
were  chosen  such  that  the  core  of  EUC-like  flow  appeared  at  the  same 
level  (~220m)  as  in  the  Newtonian  cooling  case.  (Values  of  A-K-lcm2 
sec-1  resulted  In  high  wavenumber  vertical  structure  near  the  surface 
that  bore  no  resemblance  to  observed  flow  features).  Appropriate  scales 
are  D=43m,  L=97km,  l/*llcm  sec"1,  W»4.7  x  10-3cm  sec-1,  P-2.3  x  102  dyne 
cm-2  and  T-0.22  °C;  nondimensional  zonal  wavenumber  (£q)  and  meridional 
decay  scale  (y)  are  3.0  x  10"2  and  9.3  x  10"3  respectively. 

For  the  above  winds  and  parameter  values,  bottom  trapped  solutions 

are  much  less  Important  than  surface  trapped  and  are  therefore  neglected. 

In  addition,  uzz*vzz*e2*0  at  z“0  can  &e  satisfied  by  setting  anj-an  j+g, 

j*l,2,3.  Thus  with  8  given  by  (28),  A  Is  trivially  solved  for  the  re- 

(xl 

duced  set  of  boundary  conditions  "UZZZ*T  »  VZ22“W-®  at  The  series 
expansions  for  q  were  truncated  at  n»N  -57. 


b.  Flow  Description 

Figure  23  shows  depth  profiles  on  the  equator  at  the  maximum  In 
westward  winds.  As  In  the  Newtonian  cooling  case,  one  finds  an  SEC,  EUC 
and  an  EIC  which  Is  stronger  relative  to  the  shallow  currents.  Above 
the  EUC  core  and  at  the  surface,  upwelled  water  Is  cooler  while  down- 
welling  In  and  below  the  core  Is  associated  with  warmer  water.  An  east¬ 
ward  baroclinlc  zonal  pressure  force  drives  the  EUC  and  reverses  at  depth 
to  drive  the  EIC.  Even  though  the  dlffuslvitles  are  stronger  in  this 
model,  the  vertical  wavelengths  in  Figure  23  are  decidedly  smaller  than 
In  Figure  5.  This  is  consistent  with  the  generally  smaller  depth  scales 
predicted  by  (60)  and  (61b)  relative  to  their  counterparts  (18)  and  (22b). 
On  the  other  hand,  stronger  dlffusivity  leads  to  weaker  velocity  and 
temperature  scales. 

Figure  24  Is  a  meridional  section  of  zonal  velocity  on  the  northern 
half  plane  (0<y<4)  down  to  1500m.  The  width  of  the  current  system  Is 
twice  that  shown  in  Figure  8  being  closer  now  to  the  observed  width  of 
150km.  A  tongue  of  geostrophic  eastward  momentum  sloping  downward  toward 
the  poles  emanates  from  the  EUC  much  like  the  SSCC  in  Figure  1.  A  simi¬ 
lar  tongue  of  westward  momentum  emanates  from  the  EIC. 

Meridional  circulation  in  the  upper  500m  resembles  that  shown  in 
Figure  7b  though  somewhat  reduced  In  amplitude.  A  meridional  section  of 
temperature  on  the  same  plane  (Figure  25)  shows  cooler  water  at  the  sur¬ 
face  in  the  equatorial  boundary  layer  and  doming  (troughing)  above  (below) 
the  EUC  core.  These  effects  are  in  qualitative  agreement  with  the  obser¬ 
vations  though  of  weaker  Intensity.  Moreover,  the  temperature  field  Is 
statically  stable  everywhere  in  contrast  to  that  in  Figure  22. 


Figure  23.  Depth  profiles  of  perturbation  variables  on  the  equator  for  the  biharmonlc 


Latitude  (deg) 


o  1.0  2.0  3.0 


Figure  24.  Meridional  section  of  zonal  velocity  for  the  bi 
harmonic  friction  model.  Contour  interval  is  10  cm  sec-1; 
regions  of  westward  flow  are  stippled. 
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Figure  25.  Meridional  section  of  total  temperature  for  the  biharmonic 
friction  model.  Contour  interval  is  0.33  °C.  Superimposed  are 
isotachs  of  eastward  velocity  in  the  EUC  (dashed  lines). 
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c.  Dynamics 

As  stated  earlier,  the  choice  of  boundary  conditions  (54)  guaran¬ 
tees  that  the  model  ocean  is  in  Sverdrup  balance  at  all  latitudes.  The 
flow  can  then  be  examined  in  terms  of  a  shallow  Ekman  component  and  a- 
deep  geostrophic  component.  The  balance  of  forces  in  the  Ekman  layer  is 
identical  to  that  shown  in  Figure  11  except  the  diffusive  terms  are  now 
biharmonic.  Similarly,  the  balance  for  deeper  motion  Is  given  by  Figure 
12  indicating  that  the  Kelvin  mode  is  of  central  Importance  in  the 
equatorial  boundary  layer.  As  In  Figure  14,  winds  project  onto  the 
Ekman  layer  at  mid-latitudes  and  directly  onto  geostrophic  zonal  flow  at 
low  latitudes.  The  essential  purpose  of  the  equatorial  Ekman  layer  Is 
again  to  eliminate  surface  vertical  velocities  associated  with  the  deeper 
flow. 

The  similarity  between  the  dynamics  of  this  and  the  Newtonian 

cooling  model  Is  not  surprising  since  I  have  simply  Invoked  a  different 

parameterization  for  turbulent  dissipation.  This  has  not  altered  the 

substance  of  the  force  balance,  but  merely  changed  certain  details  of 

the  flow  patterns.  For  example  in  the  northern  hemisphere,  surface 

velocity  is  68°  to  the  right  of  the  wind  in  the  mid-latitude  Ekman  layer 

as  opposed  to  45°.  The  Ekman  layer  Is  also  deeper:  De~0(y)^  as  opposed 

to  D  On  the  other  hand,  Ekman  transports  are  still  90°  to  the 

e  y 

right  of  the  wind  since  they  are  not  sensitive  to  the  form  of  the 
diffusion  operators. 

3.  Thermal  Forcing 

The  present  model  allows  for  either  perturbation  heat  flux  forc¬ 
ing  or  specification  of  surface  temperature.  I  will  briefly  discuss  so- 


lutions  for  the  latter  In  the  absence  of  wind  stress.  The  scales  D  and 
L  In  (4)  are  unchanged  but 


where  T  Is  the  temperature  scale.  The  form  of  the  temperature  boundary 
condition  Is  Identical  to  stress  condition  (27)  and  parameter  values 
are  A*Ks40cm2  sec-1,  bz310"1*  °C  cm"1.  Thus  for  T«1.0  °C,  (J»49cm  sec"1, 
W=2.2x10“2  cm  sec"1  and  P=4.3xl03  dyne  cm"2.  Series  expansions  for  q 
were  truncated  at  n»N=17  and  solved  by  the  method  described  In  Section 
C2a. 

Figure  26  shows  profiles  at  the  equator  and  at  the  longitude  of 
maximum  zonal  surface  temperature  gradient  (x=37r/4).  Colder  water  lies 
to  the  west,  warmer  to  the  east.  This  temperature  gradient  sets  up 
barocllnlc  pressure  forces  that  drive  zonal  currents  much  like  those  In 
the  wind  driven  cases.  Furthermore  temperature  and  vertical  velocity 
bear  the  same  relationship  to  these  currents  e.g.  upwelllng  (downwelllng) 
above  (below)  the  eastward  velocity  maximum  correlated  colder  (warmer) 
water.  The  width  of  the  current  system  is  again  0(L). 

The  ocean  Is  In  Sverdrup  balance,  though  all  depth  Integrated 
quantities  In  (29)  are  zero  because  t  ■  0.  The  Kelvin  mode  dominates 
the  equatorial  ocean  as  In  the  case  of  zonal  winds.  Ekman  layers 
(Figure  27)  are  needed  because  the  deep  flow  cannot  satisfy  all  the 
surface  boundary  conditions  alone.  In  the  absence  of  winds,  the 
strength  of  Ekman  layer  Is  proportional  to  the  strength  of  the  deep 
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Figure  27.  Balance  of  forces  at  the  sea  surface  for  surface  tem¬ 
perature  forcing  computed  from  just  the  Ekman  modes  In  the  (a) 
zonal  and  (b)  meridional  direction.  Forces  are  In  nondlmensional 
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flow.  This  leads  to  stronger  equatorial  rather  than  mid-latitude 
Ekman  layers  In  contrast  to  the  wind  driven  cases  (cf.  Figure  11). 

These  results  together  with  Phllander's  (1973a)  more  complicated 
equatorial  thermocllne  calculations  Indicate  that  zonal  surface  tempera¬ 
ture  gradients  can  have  a  significant  effect  on  subsurface  flow  structure. 
Models  that  attempt  to  accurately  describe  the  equatorial  current  system 
must  therefore  account  for  observed  surface  temperature  variations.  This 
will  require  a  better  understanding  of  equatorial  mixed  layer  dynamics 
and  thermodynamics  than  we  have  at  present. 

4.  Discussion 

I  have  examined  the  physics  of  a  linearized,  stratified  equatorial 
ocean  In  both  Newtonian  cooling  and  biharmonic  friction  models.  In 
either  case,  easterly  winds  excite  flows  that  resemble  the  SEC,  EUC, 

EIC  and  SSCC.  The  models  also  predict  equatorial  upwelllng  and  thermo¬ 
cllne  weakening  In  qualitative  agreement  with  observations.  The  bihar¬ 
monic  friction  model  gives  more  realistic  length,  depth  and  velocity 
scales  for  a  range  of  parameters;  it  also  provides  a  better  representa¬ 
tion  of  the  temperature  field  and  allows  surface  thermal  forcing.  The 
dynamics  of  both  models  are  identical.  In  particular,  well  behaved 
Ekman  layers  are  present  at  all  latitudes  and  the  equatorial  boundary 
layer  is  characterized  by  a  vertically  diffusive,  geostrophlc  Kelvin 
balance. 

These  models  are  distinguished  from  McCreary's  (1980)  linearized, 
stratified,  vertically  diffusive  model  by  their  ability  to  separate  and 
follow  the  development  of  Ekman  and  geostrophlc  layers  as  they  approach 
the  equator.  McCreary  employs  diffusion  operators  which,  unlike  the 
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Newtonian  cooling  and  biharmonic  friction  formulations,  lead  to  a 
Sturm-Llouvllle  eigenvalue  problem  for  the  vertical  dependence.  The 
elgenmodes  either  Individually  or  In  groups  cannot  be  Identified  with 
Ekman  or  geost'rophlc  dynamics  as  was  done  with  the  exponential  modes  4 
In  Section  Cl.  Moreover,  the  biharmonic  friction  model  allows  surface 
thermal  forcing  which  was  not  allowed  In  McCreary's  model.  However 
these  differences,  resulting  from  different  turbulence  parameterlza- 
tlons  do  not  obscure  basic  similarities  between  the  two  models  which 
have  been  coimented  on  throughout  the  text.  Thus,  this  and  McCreary's 
(1980)  work  should  be  viewed  as  complementary  approaches  to  under* 
standing  steady  equatorial  dynamics. 

To  stimulate  discussion,  extensive  comparison  has  been  made 
between  model  results  and  observations.  However  the  purpose  of  this 
study  has  not  been  to  simulate  ocean  flows.  Rather  It  Is  hoped  that 
the  results  of  this  work  will  clarify  Important  processes  In  more  com¬ 
plicated  numerical  models  and  also  motivate  further  modeling  of  the 
equatorial  ocean. 


CHAPTER  III 

EQUATORIAL  KELVIN  AND  I NERTI O-GRAVITY  WAVES 
IN  ZONAL  SHEAR  FLOW 

Theories  of  the  equatorial  ocean  circulation  have  evolved  along 
two  separate  and  distinct  lines.  In  Chapter  II  I  discussed  the  older 
of  these,  viz.  frictional  models  of  the  mean  flow  which  date  from  the 
mid  1950's.  While  they  contribute  to  our  understanding  of  steady 
equatorial  dynamics,  an  obvious  limitation  of  such  models  is  their  in¬ 
ability  to  account  for  the  wide  variety  of  temporal  fluctuations  found 
In  modern  oceanic  measurements.  To  overcome  this  limitation,  linear, 
time  dependent,  Invlscld  equatorial  wave  theory  was  developed  in  the 
early  1960's.  Wave  theory  has  since  provided  satisfying  descriptions 
of  heretofore  puzzling  oceanic  phenomena  such  as  the  Initiation  of  an 
El  Nino  (McCreary,  1977),  Somali  Current  dynamics  (Llghthill,  1969; 
Anderson  and  Rowlands,  1976)  and  sea  level  variability  (Wunsch  and  Gill, 
1976).  However  the  basic  theory  Is  formulated  for  small  amplitude 
perturbations  In  an  ocean  with  no  mean  currents.  Considering  the  ob¬ 
served  strength  (-100  cm  sec-1)  of  the  EUC  and  SEC,  it  Is  likely  that 
waves  and  steady  currents  interact  In  manner  that  neither  wave  nor  current 
theory  alone  can  predict.  This  chapter  Is  one  of  the  first  attempts  to 
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Include  In  an  equatorial  wave  model  mean  flows  like  those  observed  In 
the  tropical  ocean. 
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ABSTRACT 

The  interaction  of  fcottrophic  zonal  mean  currents  in  an  equatorial  ocean  with  free,  neutrally  stable, 
internal  Kelvin  and  inert io- gravity  waves  is  investigated  using  a  two-layer,  reduced-gravity  model.  Solu¬ 
tions  in  the  inertio-gravity  range  are  obtained  by  a  simple  numerical  integration  scheme  that  allows 
several  dilferem  background  flows  to  he  tested.  It  is  found  that,  due  to  the  interaction,  the  amplitudes 
and  latitudinal  distributions  of  wave  zonal  velocity  are  substantially  altered  from  those  which  would 
occur  in  the  absence  of  mean  flow.  Meanders  of  currents  similar  to  those  observed  during  GATE  may 
be  interpreted  as  advcctions  of  mean  currents  by  wave  meridional  velocity.  On  the  other  hand,  wave 
pressure  (sea  level)  and  meridional  velocities  are  not  greatly  affected  by  the  mean  flow.  These  results 
may  be  of  importance  in  attempting  to  At  equatorial  wave  theories  to  observations  of  zonal  current, 
on  the  one  hand,  and  of  sea  level  fluctuations  on  the  other. 

Kelvin  waves  are  treated  using  a  perturbation  expansion  based  on  the  small  ratio  of  mean  current 
speed  to  wave  phase  speed.  The  shear  flow  alters  low-frequency  Kelvin  waves  only  slightly,  intro¬ 
ducing  a  small  meridional  velocity  and  a  Doppler  shift  which  could  atfect  the  speed  of  baroclinic  ad¬ 
justments  in  the  tropics.  At  higher  frequencies,  the  Kelvin  wave  becomes  more  like  an  inertio-gravity 
wave  and  may  appear  in  velocity  records  as  varicose  meanders  of  the  background  current. 


I.  Introduction 

The  theory  of  time-dependent  motion  in  the 
tropics  is  distinguished  from  midlatitude  theories  by 
the  fact  that  the  Coriolis  parameter  changes  sign  at 
the  equator,  creating  a  waveguide  in  which  there 
can  be  no  net  meridional  propagation  of  energy. 
Within  the  waveguide,  two  types  of  motion  unique 
to  the  tropics  are  found  (Fig.  i).  First,  internal  east¬ 
ward-traveling  Kelvin  waves  can  exist  in  the  ab¬ 
sence  of  physical  boundaries.  Because  their  group 
velocity  is  the  fastest  of  any  open  ocean  internal 
wave,  they  afford  a  most  efficient  mechanism  for 
the  zonal  transfer  of  energy  and  momentum,  sub¬ 
stantially  reducing  baroclinic  adjustment  times  at 
low  latitudes.  For  example,  both  McCreary  ( 1976) 
and  Hurlburt  et  at.  (197b)  found  in  simulations  of 
El  Nino  that  when  the  easterlies  relax  in  the  tropics, 
internal  Kelvin  waves  excited  at  the  western  bound¬ 
ary  can  cross  a  basin  several  thousand  kilometers 
wide  in  1-2  months.  Various  authors  have  shown 
that  wind-generated  Kelvin  waves  in  the  western 
Atlantic  are  capable  of  producing  the  seasonal  up- 
welling  signals  observed  along  the  coast  of  Guinea, 
at  the  opposite  side  of  the  ocean  ( Moore  et  al. .  1978: 
Adamec  and  O'Brien.  1978).  These  and  other  studies 
suggest  that  Kelvin  waves  should  be  a  prominent 
feature  of  the  equatorial  circulation,  although  none 
have  been  detected  in  the  ocean  (Wunsch.  1977a). 

The  second  class  of  motions  shown  in  Fig.  I  con- 
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sists  of  long  inertio-gravity  waves  which  for  low- 
order  vertical  modes  typically  have  periods  between 
-2  days  and  -2  weeks  and  wavelengths  >2500  km. 
This  group  of  intermediate  frequency  disturbances 
is  of  geophysical  interest  because  it  includes  waves 
which  by  virtue  of  their  small  zonal  group  velocity 
are  preferentially  excited  by  random  local  forcing. 
It  also  includes  the  mixed  Rossby-gravity  or  Yanai 
wave.  Matsuno  ( 1966)  argued  that  because  of  the  ex¬ 
istence  of  such  free  modes,  one  should  in  principle 
expect  more  energy  at  periods  between  a  few  days 
and  a  few  weeks  in  the  open  equatorial  ocean  than 
at  midlatitudes.  Many  investigators  have  found  this 
energy  (Wunsch.  l9T7a).  but  to  identify  it  with  par¬ 
ticular  waves  is  difficult  because  most  studies  have 
lacked  the  necessary  spatial  and/or  temporal  resolu¬ 
tion.  One  notable  exception  is  the  Wunsch  and  Gill 
1 1976)  study  of  Pacific  sea  level  variability.  Using 
a  data  base  consisting  of  long  sea  level  and  weather 
records  from  widely  scattered  island  stations,  they 
concluded  that  first  baroclinic  mode  equatorial  in¬ 
ertio-gravity  waves  excited  by  local  winds  were  re¬ 
sponsible  for  the  significant  energy  found  in  the  •*- 
day  and  5-day  bands.  Their  calculations  necessarily 
emphasized  spectral  resolution,  leaving  consider¬ 
able  amplitude  uncertainty  in  their  results. 

Aside  from  the  lack  of  data,  a  fundamental  prob- 
•em  in  interpreting  oceanic  phenomena  in  terms  of 
equatorial  waves  lies  in  the  limitations  of  the  linear 
theory  itself.  Specifically,  little  is  known  about  how 
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these  waves  interact  with  highly  sheared  mean  cur¬ 
rents  such  as  those  encountered  in  the  tropical 
oceans  since  most  of  the  basic  theory  of  time-de¬ 
pendent  equatorial  motion  (e.g..  Blandford.  1966; 
Matsu  no,  1966:  Moore.  1968:  Moore  and  Philander, 
1977)  considers  perturbations  linearized  about  a 
state  of  rest.  Existing  theories  of  wave-mean  shear 
interaction  in  the  atmosphere  (Boyd,  1978)  cannot 
be  applied  directly  to  the  ocean  since  they  depend 
on  perturbation  analyses  in  which  the  small  param¬ 
eter  is  the  ratio  of  vertical  or  meridional  wave  scales 
to  corresponding  mean  shear  scales;  in  the  ocean, 
this  ratio  is  CXI)  or  larger. 

The  purpose  of  this  paper  is  to  investigate  the 
Kelvin  wave  and  the  long,  intermediate-frequency 
inertio-giavity  waves  in  the  presence  of  strong,  merid¬ 
ionally  sheared  geostrophic  jets  appropriate  to  the 
equatorial  ocean.  We  treat  these  disturbances  in  a 
two-layer  model  in  which  zonal  background  flows 
are  confined  to  the  upper  layer  in  which  the  wave 
scales  ate  comparable  to  those  of  the  first  baroclinic 
•node  in  a  realistically  stratified  ocean.  The  analysis 
for  this  system  is  especially  simple  because  the 
phase  speeds  are  much  greater  than  typical  mean 
current  speeds  so  that  there  are  no  critical  layers 
and  no  growing  instabilities.  Slower  waves,  i.c., 
those  associated  with  higher  vertical  modes,  could 
be  'instable,  but  we  do  not  consider  these.  Further¬ 
more.  because  we  are  not  attempting  to  simulate  a 
particular  set  of  observations  but  rather  to  isolate 
the  structural  and  dispersion  modifications  due  to 
mean  currents  in  a  few  of  the  gravest  horizontal 
modes,  no  forcing  is  imposed  on  the  system. 

In  related  work  Philander  ( 1976:  1978a)  adapted 
the  discrete  layer  model  to  a  discussion  of  critical 
layers  in  highly  sheared  zonal  flows  along  the 
equator.  He  found  that  baroclinic  instabilities  are 
suppressed  in  the  tropics  and  that  westward  flows 
are  more  barotropically  unstable  than  eastward.  For 
typical  values  of  oceanic  parameters,  the  fastest 
growing  disturbances  propagated  westward  with 
periods  in  excess  of  three  weeks  and  wavelengths 
shorter  than  2000  km.  He  also  found  solutions  cor¬ 
responding  to  longer,  neutrally  stable  waves,  but 
these  were  not  discussed  in  detail. 

Hallock  (1977)  used  a  similar  three-layer  numeri¬ 
cal  model  with  CX40  cm  s~')  mean  currents  in  the 
upper  two  layers  to  simulate  meanders  of  the  under¬ 
current  and  South  Equatorial  Current  observed  dur¬ 
ing  Phase  II  of  the  GARP  Atlantic  Tropical  Experi¬ 
ment  (GATE)  (Doing  el  at.,  1975).  He  impulsively 
stressed  the  surface  of  the  ocean  at  a  fixed  (2400 
km)  wavelength  with  a  northerly  wind  stress 
anomaly,  thereby  exciting  a  neutrally  stable.  Yanai- 
Uke  disturbance  that  advected  the  mean  currents 
back  and  forth  across  the  equator  with  the  observed 
16-day  period  and  zonal  phase  propagation.  Higher 
frequency  gravity  waves  appeared  as  noise  at  the 
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same  wavelength:  Kelvin  waves  were  not  excited 
due  to  the  symmetry  of  the  forcing. 

Most  recently.  Philander  (1979)  has  considered 
neutrally  stable  Kelvin  and  inertio-gravity  waves  in 
a  two-layer  numerical  model  with  eastward  shear 
flow  in  the  upper  layer.  For  typical  values  of  oceanic 
parameters,  he  finds  that  Kelvin  waves  are  Doppler- 
shifted  to  higher  frequency  and  have  meridional  ve¬ 
locities  which  increase  with  increasing  frequency. 
This  meridional  velocity  does  not  affect  the  geo¬ 
strophic  balance  of  the  zonal  currents  except  at  pe¬ 
riods  <8  days,  beyond  which  the  waves  become 
more  like  inertio-gravity  oscillations.  He  also  found 
that  the  dispersion  of  inertio-gravity  waves  is  little 
changed  by  the  presence  of  eastward  currents,  but 
details  of  the  wave  structures  were  not  discussed. 

There  is  overlap  between  the  Philander  (1979) 
study  and  the  one  presented  here.  Both  consider 
inertio-gravity  waves  using  numerical  techniques, 
but  Philander's  analysis  covers  a  wider  range  of 
wavenumbers.  On  the  other  hand,  by  restricting  our 
attention  to  a  limited  region  near  zero  wavenumber, 
we  have  simplified  the  numerical  integration  and 
are  able  to  treat  in  more  detail  the  behavior  of  these 
waves  in  several  different  flow  regimes.  Regarding 
the  Kelvin  wave,  we  exploit  the  fact  that  typical 
mean  flow  speeds  are  much  less  than  the  phase 
speeds  of  low-order  vertical  modes  in  realistically 
stratified  oceans.  Using  the  small  ratio  of  these 
speeds  in  a  perturbation  expansion,  we  obtain  re¬ 
sults  similar  to  Philander's  1 1979)  numerical  results 
but  in  which  both  eastward  and  westward  mean 
flows  are  considered.  Philander  in  the  same  paper 
also  considered  the  low-frequency  Rossby  waves  in 
eastward  flow,  a  subject  not  treated  here. 

In  Section  2.  we  develop  the  model  equations  for 
a  two-layer  ocean  with  geostrophic  mean  currents 
in  the  upper  layer.  Section  3  describes  the  numeri¬ 
cal  solution  of  these  for  intermediate  frequency  in¬ 
ertio-gravity  waves  at  and  near  zero  wavenumber. 
The  most  striking  result  of  this  section  is  the  en¬ 
hanced  variability  of  zonal  velocity  fluctuations  in 
mean  shear  zones.  In  Section  4.  we  describe  the 
perturbation  expansion  for  the  Kelvin  wave  which 
yields  Doppler-shifted  frequencies  and  nonzero 
meridional  velocities  proportional  to  a  variety  of 
flow  parameters.  Section  5  presents  a  discussion  of 
the  results  of  Sections  3  and  4. 

2.  Model  ocean 

We  consider  a  two-layer,  in  viscid.  Boussinesq 
ocean  on  an  equatorial  /3-plane  with  no  lateral 
boundaries.  We  use  a  caret  (omitted)  to  denote  the 
dimensional  (nondimensional)  form  of  a  variable. 
In  the  upper  layer  of  constant  density  flows  a  steady 
zonal  current  0  varying  in  latitude  only  and  geo- 
strophically  balanced  by  6,  the  thickness  of  the 
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layer.  The  lower  layer  of  density  p,  >  p»  is  infinitely 
deep  and  therefore  motionless.  Using  nondimen- 
sionalizations  explained  in  (3)-(4)  and  linearizing 
about  this  background  flow,  the  vertically  integrated 
equations  governing  perturbations  of  the  form 
etw'-w'aie 

-iau  +  R (ikllu  +  vU,)  -  yv  +  ikp  -  0,  (la) 


r-^-T -17BF'  (4) 

which  measures  the  ratio  of  background  flow  speed 
to  the  phase  speed  of  internal  gravity  waves. 

By  setting  R  *  0.  00  *  I  and  A2  =  g'D,  we  ob¬ 
tain  the  well-known  equations  for  perturbations 
about  a  state  of  rest: 


-itrv  ♦  iiR Uv  +  yu  *  p»  “  0,  (lb) 

A* 

— -r  ( ikKU  -  ia)p  +  ikDu  *■  (Du),  »  0,  (Ic) 

go 

where  r  is  positive  eastward,  y  positive  northward 
ami  the  reduced  gravity  g'  «  tg  «  g{ !  -  pjp,)- 
The  set  (u.r.p)  gives  the  meridional  dependence 
of  zonal  velocity,  meridional  velocity  and  pressure 
fluctuation,  respectively.  The  background  current 
0  and  layer  thickness  D  are  given  by 

D  -  flD(y)  - f)[D,  +  RD,(.v)],  D,(-*)-0.  (2a) 

|-7g-jyC/(.v)  -  -D„ty),  (2b) 

where  D  is  the  scale  of  upper  layer  thickness  and 
D,  is  an  0(1)  constant  and  accounts  for  small 
differences  in  avetage  layer  thickness  between  the 
several  background  flows  to  be  studied.  The  non- 
dimcnsionalizations  and  parameter  definitions  used 
to  obtain  (l)-(2)  are 

(i.y)  -  lu,y),  H  •  k/L,  L  -  (/3/A)-*'* 
i  iT.  a  -  trlT,  T  «  (/8A)'1- 
0  m  OU.  O  »  max|(/| 

(«,£,(p4Ar'p)  =»  -&(  M.c.p) 

The  parameter  A  appears  in  ( I )— (3)  as  an  artifact 
of  the  nondimensionalization.  ensuring  an  inverse 
relation  between  length  and  lime  scales.  In  fact, 
it  is  the  parameter  which  would  have  arisen  as  the 
separation  constant  or  eigenvalue  of  the  vertical 
problem  had  we  not  proceeded  ab  initio  to  the  two- 
layer  limit  but  retained  a  continuous  system  in  the  ver¬ 
tical.  It  becomes  important  to  recognize  the  eigenvalue 
nature  of  A  in  Section  4.  There  we  treat  the  Kelvin 
wave  by  perturbation  expansions,  and  as  is  usual 
in  such  problems  it  will  be  found  necessary  to  allow 
perturbations  to  A  in  order  to  avoid  spurious  reso¬ 
nances  at  first  order  driven  by  zeroth-order  solu¬ 
tions.  In  this  section  and  Section  3,  where  direct 
numerical  calculations  are  used,  we  simply  set  A* 
•  g‘t>  throughout. 

The  parameter  V  is  a  scale  for  the  perturbations 
and  is  essentially  arbitrary,  provided  >/**£/.  so  that 
the  basic  linearization  is  valid.  The  remaining  pa¬ 
rameter  in  (l)-(3)  is  the  Rossby  number 


vn  +  (a1  -  k*  -  kto  -  y2)v  «  0,  (5a) 

u  -  i(«r*  -  k*)~Hayv  -  kv,),  (5b) 

p  -  i(«r*  -  k')-<(kyv  -  crt>,).  (5c) 

Solutions  to  (5),  decaying  as  |  y  |  -»  ®,  exist  for 

a*  -k1  -  klo  «  2n  t  1,  n  -  1,2 .  (6a) 

with 

o  -  *.(y)  «  <ir‘'22*n'.)-'’  H.(y)e-It,  (7) 


where  W.(v)  are  Hermite  polynomials.  The  high- 
frequency  limit  of  (6a)  yields  inertio-gravity  waves, 
whereas  the  low-frcquencv  limit  yields  Rossby 
waves.  Turning  latitudes,  where  ihe  solutions 
change  from  oscillatory  to  exponential  behavior,  are 
located  at  y  *  r(2«  +  l)‘J. 

For  n  =  0  the  only  root  of  ;6ai  which  gives  a 
solution  with  bounded  u  and  p  fields  at  >  =  -*■  is 

k  •  <r - -  .  (6b) 

a 

This  is  the  mixed  Rossby-gravitv  or  Yanai  wave  of 
Fig.  I  which  for  it  ►  I  is  gravity-like  and  o  <*  1  . 
is  Rossby-like. 

One  more  solution  to  (I)  when  R  =  0  which  is 
not  contained  in  1 5)  is  the  Kelvin  wave  (designated 
n  *  -I  in  Fig.  I),  given  by 

v  -  0.  (8a) 

U  m  p  =  tfl,.  (gb) 

<r  -  k.  (8c) 

The  remainder  of  this  paper  deals  with  solutions 
to  (I)  when  R  *  0.  We  fix  the  parameters  t  «  4 
x  I0'2  and  D  =  200  m  so  that  the  Rossby  radius 
(Z.  »  360  km),  time  scale  (T  »  1 .46  days)  and  eigen¬ 
value  (A  -  280  cm  s_l)  correspond  closely  to  those 
of  the  first  barociinic  mode  in  a  realistically  strati¬ 
fied  ocean  (Moore  and  Philander.  1977).  We  note 
that  the  choice  of  D  is  a  compromise.  A  200  m  upper 
layer  coincides  approximately  with  the  depth  range 
of  the  strongest  mean  currents.  It  also  corresponds 
to  the  layer  of  greatest  horizontal  velocity  ampli¬ 
tudes  in  the  first  barociinic  mode  computed  for  a 
realistically  stratified  resting  ocean  (Wunsch.  1975: 
Philander.  1978b).  However,  the  node  of  the  first 
barociinic  mode  thus  computed  lies  at  about  1200  m 
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Fk>.  I.  Nondifflcnsionalind  dispersion  diagram  for  equatorial  /3-plane  waves. 
Stippled  region  is  discussed  in  teat.  (After  Wunsch  and  Gilt.  1976.) 


not  200  m  and  to  this  extent  the  two-layer  model 
differs  from  the  true  vertical  structure. 

For  the  inenio-gravity  waves  treated  numerically 
in  the  following  section,  the  above  choice  of  param¬ 
eters  requires  that  each  of  the  different  background 
states  examined  satisfies 

t>  «  (l/2y.)  I*"  bdy  -  200  m,  (9) 

with  D,  in  (2)  determined  from  the  integration  con¬ 
stant.  The  range  (-y„y0)  is  the  range  of  significant 
wave  amplitude,  taken  as  ±13*  of  latitude.  Since 
the  Kelvin  wave  discussed  in  Section  4  is  solved 
using  perturbation  expansions,  we  can  standardize 
those  results  to  the  specified  (c.f))  by  defining  a 
convenient  zeroth-order  mean  mixed  layer  depth. 

3.  lnertio-gravity  waves 
a.  Method  of  solution 

We  set  A*  »  g't>  and  combine  (1)  to  gel  expres¬ 
sions  for  ■>,  u,  p: 

Do„  -  iyRUv,  -  [**£>  -  o*  +  y*  +  RU 
+  *£>*"(1  -  -  (H)-' 

x  (2okRU,  -  k*yRU)[kDv{y  -  RU.) 

-  UDo,  -  yvRU )]  -  0,  (10a) 
«  -  i{k'D  -  d*)-'[*(0i-,  -  y»R If) 

-  <7tK>  -  Rt/,)],  (10b) 
p  -  i{k*D  -  o’J-'totOv,  -  yvRU) 

-  kDvky  -  Rt/,)],  (10c) 
&mV-  kRUly). 


We  restrict  our  attention  to  solutions  of  <  10a)  for 
the  stippled  region  in  Fig.  I  bounded  by  |  k  J  <  1.0 
and 0.7  <  a  <  3.0  (which  dimensionally  corresponds 
to  wavelengths  >2500  km  and  periods  between  3 
and  13  days)  in  order  to  avoid  singularilies  ti.e.. 
critical  layers)  atir-ARC/  =  0  and  apparent  singu¬ 
larities  at  Dk *  -  Ur  -  k  R  t/)’  =  0.  This  not  only  re¬ 
duces  the  amount  of  computation  needed  to  obtain 
v  but  guarantees  all  the  wave  solutions  will  be  neu¬ 
trally  stable. 

Eq.  (10a)  was  integrated  using  a  simple  fourth- 
order  Runge-Kutta  routine  by  shooting  for  the  eigen¬ 
value  a  given  a  choice  of  k.  We  used  boundary 
conditions  on  i>  similar  to  those  used  by  Philander 
(1976.  1978a.  1978b)  and  Hallock  ( 1977).  viz.,  v  -  0 
at  y  »  ♦<*>.,  v,c,  continuous  at  y  -  0.  For  sym¬ 
metric  U.D  there  is  the  added  condition  at  y  «  0 
that  v  or  v,  be  zero  depending  on  the  mode  under 
consideration.  Integration  proceeded  from  the  pole¬ 
ward  boundary  to  eliminate  solutions  to  (10a)  that 
grow  as  |  y  |  -*  ».  To  test  the  accuracy  of  this  pro¬ 
cedure,  numerical  solutions  to  (10)  for  R  -  0  and 
step  size  A  =*  0.02  were  compared  to  the  known 
analytic  solutions  of  (5)  given  by  (7).  Table  1  shows 
that  except  near  the  boundaries  y  »  ±6,  one  can 
expect  four  decimal  place  accuracy  in  both  eigen¬ 
values  and  eigenfunctions.  This  accuracy  is  assumed 
when  (10)  is  integrated  for  R  +  0. 

b.  Results 

Eq.  (1)  was  integrated  for  four  geostrophic  back¬ 
ground  states,  three  of  which  are  given  by 

RU  -  O  Se-’**. 

D  -  0.99  +  0.083e-»-,  R  -  0.5,  (11a) 
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Tails  I.  Accuracy  of  numerical  solutions  (subscript  A)  compared  to  analytic  solutions  for  t  -  0,  step  size  A  »  0.02. 


y 

V 

«A  « 

Pb 

P 

0. 

-0.07259 

-0.07250 

0.  0. 

0.93290  x  I0-* 

0. 

1.0 

0.52922 

0.52925 

0.23660  0.23667 

-0.70993 

-0.71006 

2.0 

0.02665 

0.52664 

0.73940  0.73930 

0.31675 

0.31607 

3.0 

0.16479 

0.16479 

0.22110  0.22100 

0.16909 

0.16906 

4.0 

0.90730  x  10-* 

0.90743  x  10-' 

0.16234  x  IO-‘  0.16232  x  10  ' 

0.14147  x  |0-‘ 

0.14130  x  10" 

5.0 

0.15934  x  10-* 

0.15933  x  I0-* 

0.35631  x  I0-*  0.35627  x  I0-* 

0.32765  x  |0-> 

0.32720  x  10-* 

5.9 

0.96765  x  10-* 

0.16533  x  I0-* 

0.25532  x  I0-*  0.43624  x  10-* 

0.57629  x  |Q~» 

0.41001  x  t0-> 

«»■  -  4.9997  a*  -  5.0000 

Rtf  -  - 0.5e -*»*, 

D  -  1.01  -  0.083c-J‘',f  R  -  0.5,  (11b) 


tion  3b2.  The  dispersion  of  inertio-gravity  waves 
in  ail  four  mean  fields  is  analyzed  in  Section  3b3. 


R U  »  -ye-", 

D  -  0.56  -  0.5  ^ye -  j'  e"Vf  j  , 
R-0.43.  (11c) 


1)  Meridional  structure:  k  «  0 
When  k  «  0,  Eq.  ( 10)  becomes 
Dv„  -  2yRl/v,  +  (o*  -  y‘  -  RU)v  »  0,  (12a) 
u  «  i<r~'(y  -  R U,)v.  (12b) 


The  flows  of  (1  la)  and  ( 1  lb)  are  symmetric  eastward 
and  westward  jets,  respectively,  with  maximum 
speeds  of  140  cm  s~‘  or.  the  equator  and  half-width 
of  approximately  1.5°.  The  flow  in  (lie)  is  an  anti¬ 
symmetric  jet  with  maximum  speed  of  120  cm  s~' 
at  ±  2.3°.  In  each  case  the  mixed-layer  depth  is  speci¬ 
fied  by  (2)  and  (9)  with  y0  =  4  (dimensional  latitude 
=  13").  These  three  background  slates  are  illustrated 
in  Fig.  2.  Note  in  parti:ular  the  very  deep  mixed 
layer  to  the  north  associated  with  the  antisymmetric 
jet  (lie). 

Also,  we  integrated  ( 10)  for  a  fourth,  less-idealized 
background  state,  derived  from  a  temperature  sec¬ 
tion  typical  of  average  conditions  in  the  upper  layers 
of  the  central  Pacific  (Fig.  3  (after  Wyrtki  ct  at., 
1977)1.  We  chose  the  interface  between  the  constant 
density  layers  of  the  model  to  be  the  20°C  isotherm 
which  is  centrally  located  in  the  equatorial  pycno- 
cline;  a  constant  of  40  m  was  added  to  the  mean 
depth  as  determined  by  the  isotherm  so  as  to  satisfy 
b  -  200  m.  Assuming  that  gcostrophy  holds  to 
within  ±0.5°  of  the  equator,  we  calculated  a  current 
system  in  balance  with  this  interface  by  cubic  spline 
interpolation  of  ( -D,Jy ).  The  latitudinal  distribu¬ 
tion  (Fig.  3,  dashed  line)  and  mass  transport  of  the 
interpolated  zonal  currents  are  in  agreement  with 
observations.  For  example,  the  model  undercurrent 
transport  of  -30  x  10“  m3  s*1  compares  favorably 
with  Knauss'  (1966)  values.  The  Rossby  number  for 
this  case  is  -0.25.  Choosing  other  isotherms  in  the 
thermocline  to  define  the  model  interface  does  not 
significantly  change  this  Pacific  profile. 

In  Section  3b I  we  consider  solutions  to  (10)  when 
k  «•  0  for  the  background  states  (II)  and  for  the 
Pacific  profile  of  Fig.  3.  Next  we  discuss  solutions 
to  ( 10)  for  small  but  nonzero  wavenumbers  in  Sec- 


p  «  -i<r~‘lDv),.  (12c) 

Figs.  4-6  show  solutions  to  ( 12)  in  selected  back¬ 
ground  flows  superimposed  on  solutions  corre¬ 
sponding  to  R  =  0  (dashed  lines):  the  background 
currents  themselves  are  plotted  for  reference  f  dotted 
lines).  All  solutions  have  been  arbitrarily  scaled  so 
that  the  maximum  meridional  velocity  is  1. 

Fig,  4  is  a  plot  of  solutions  to  1 12)  for  the  east¬ 
ward  jet  ( I  la).  The  effects  of  this  mean  current  on 
v  and  p  are  minimal  because  terms  involving  R  are 
relatively  small  in  ( 1 2a)  and  ( 1 2c).  On  the  other  hand, 
the  zonal  perturbation  velocity  exhibits  large  anoma¬ 
lies  in  the  region  of  strong  mean  shear,  which  from 
(12b)  is  clearly  due  to  meridional  advection  of  the 
mean  current.  The  net  result  is  a  sinuous  ( v  even) 
or  varicose  ( v  odd)  meander  of  the  background  flow. 
Similar  meanders  associated  with  the  mixed  Rossby- 


Fio.  2.  Profiles  of  upper  layer  thickness  and  associated  zonal 
current  for  three  of  the  background  flow  regimes  studied:  la) 
is  the  symmetric  eastward  jet.  tbt  the  symmetric  westward  jet 
and  tc)  the  antisymmetric  jet. 
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latitude 

Fio.  J.  Temperature  section  from  the  Central  Pacific.  The  20*  isotherm  (heavy  solid  line)  defines  the 
interface  topography  and  currant  structure  (dashed  line)  of  the  Pacific  profile  discussed  in  test. 

gravity  wave  were  the  dominant  motions  in  Hal-  teger  constants  (2n  +  I)  predicted  by  (6).  bui  are 
lock's  (1977)  wind-driven  model,  but  at  slightly  fractionally  larger  (Table  2).  Nonetheless,  solutions 
higher  wavenumber  and  lower  frequency.  Note  that  to  (10>  and  (12)  will  still  be  designated  b>  t,  -  0.  1.2. 
the  eigenvalues  a1  of  (12a)  are  no  longer  simple  in-  Not  shown  are  solutions  to  (12)  for  the  westward 


Fio.  4.  Comparisons  of  wevt  fialds  with  and  without  bucL|round  flow  tin  this  instance 
Iht  eastward  jet),  for  *  -  0.  Mendional  modes,  0.  I.  2  ere  illustrated.  Dished  Isotidl  line 
shows  the  wave  amplitude  without  (with)  background  flow.  The  background  currant  itself 
it  indicated  by  the  additional  dolled  curve. 
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Tama  2.  T he  efftctt  of  mean  tonal  flow  oo  frequency  el  tero  wavenumber.  The  percent  of  chanfe  in  frequency 

it  defined  at  IOOKo/Om)  -  I). 


Mean  flow 
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1.00 

9.17 

1 
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-1 
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3.29 
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3 

jet  (lib)  since  changing  the  sign  of  RU  essentially 
changes  the  sign  of  the  structural  and  frequency 
anomalies  observed  in  Fig.  4. 

The  antisymmetric  jet  tile)  produces  more  dra- 


RU«-»exp(-y*) 


Pro.  5.  At  in  Fi|.  4  except  with  antisymmetric  jet  at  background  flow. 


matic  alteration  of  the  eigensolutions  to  ( 12).  as 
shown  in  Fig.  3.  With  increasing  it.  perturbation 
velocities  weaken  in  the  hemisphere  of  deepest 
mixed  layer,  especially  near  the  turning  latitude 
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(equal  to  the  inertial  latitude  when  k  ■  0)  located 
tty  *  (2n  +  I)1'1.  Meridional  advection  of  the  mean 
flow  distort*  the  zonal  perturbation  velocities  even 
more,  patticulariy  that  of  the  mixed  Rossby-gravity 
wave  which  at  the  equator  shows  a  near  maximum 
instead  of  a  node.  Boyd  ( 1978)  reports  a  similar  fea¬ 
ture  associated  with  long,  atmospheric  mixed  mode 
waves  in  weakly  sheared  mean  winds  having  an 
antisymmetric  component.  Pressure  fluctuations  are 
weaker  relative  to  their  strength  when  no  mean  cur¬ 
rents  are  present,  but  their  meridional  structures 
are  otherwise  little  changed. 

Results  for  the  Pacific  profile  are  shown  in  Fig.  6. 
As  in  the  case  of  the  antisymmetric  jet,  strong  mean 
shear  right  on  the  equator  results  in  large  zonal  ve¬ 
locity  fluctuations  in  the  a  -  0  mode.  However, 
since  the  bad$g?«und  flow  between  the  turning  lati¬ 
tudes  is  largely  symmetric,  there  is  little  trend 
toward  consistently  higher  velocities  in  one  hemi¬ 
sphere  relative  to  the  other,  it  is  expected  that  such 
trends  would  appear  in  higher  meridional  modes 
whose  turning  latitudes  are  farther  from  the  equator. 
Other  effects  due  to  the  presence  of  mean  flow  in¬ 
clude  the  small-scale  variability  introduced  into  the 
latitudinal  structure  ofu.  which  corresponds  to  sim¬ 
ilar  small-scale  variability  in  the  mean  shear  pro¬ 


file.  Pressures  are  again  relatively  insensitive  to  the 
background  state. 

Since  there  are  no  critical  layers  in  the  flow  fields 
described  above,  perturbation  energy  density  aver¬ 
aged  over  one  wave  period  must  be  independent 

of  time,  i.e., 

£  -  £(y)  -  K(y)  +  V(y).  (13a) 

£(y)  -  V40(s?  +■?),  (13b) 

V(y)  -  Wp,  (13c) 

where  £  is  the  total  energy  density,  K  the  kinetic 
energy  density  and  V  the  potential  energy  density. 
The  overbar  denotes  time  average.  Furthermore, 
when  the  background  state  is  a  symmetric  function 
of  latitude,  so  must  £,  K  and  V'  be  symmetric.  This 
does  not  hold  in  the  case  of  the  antisymmetric  jet 
(lie)  or  the  Pacific  profile.  For  example,  plots  of 
K  and  V  versus  latitude  for  the  modes  n  *  0,  1.2 
when  the  mean  state  is  defined  by  (lie)  show,  as 
expected,  that  kinetic  energy  density  is  amplified 
in  the  hemisphere  of  shallow  mixed  layer  (Fig.  7). 
The  most  significant  asymmetry  appears  near  the 
turning  latitudes  except  in  the  a  -  0  mode  which 
is  complicated  by  large  zonal  perturbation  velocities 
near  the  equator.  The  dashed  curves  in  Fig.  7  that 


Pacific  profile 
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Fio.  7.  Kiiwic  ( K )  and  potential  ( V)  entity  deniitiee  aa  functiom  of  latitude,  both  with  (tolid 
Una)  and  without  (dashed  line)  borhtround  tow,  for  the  same  modes  as  Fig.  3.  The  background 
low  is  the  antisymmetric  jet. 


correspond  to  K  for  a  mean  state  of  rest  display 
none  of  these  features.  Conversely,  potential  enerxy 
density  is  similar  to  that  in  the  absence  of  mean 
currents.  (The  large  discrepancy  in  total  energy 
levels  which  appears  in  Fig.  7c  is  not  significant 
because  we  have  arbitrarily  scaled  the  amplitude 
of  the  solutions  according  to  -  I  for  each  in¬ 
dividual  wave.) 

2)  Meridional  structure:  |*|  <  1 

Solutions  to  (10)  at  small  but  non-zero  wavenum¬ 
bers  show  the  same  basic  wave-mean  flow  interac¬ 
tions  discussed  above.  Specifically,  in  regions  where 
mean  shear  is  large,  zonal  perturbation  velocities 
are  stronger;  where  the  background  fields  are  not 
symmetric  between  turning  latitudes,  kinetic  energy 
density  favors  the  hemisphere  of  shallow  mixed 
layer;  perturbation  pressures  and  meridional  veloci¬ 
ties  are  typically  less  affected  except  for  an  ampli¬ 
tude  factor.  Because  of  the  similarity  of  these  re¬ 
sults  with  zero  wavenumber  solutions,  we  do  not 
present  them  graphically. 

3)  Dispersion 

Fig.  8  compares  the  dispersion  curves  for  selected 
modes  with  and  without  background  Row.  In  Fig.  8a 
there  is  an  increase  in  eastward  and  decrease  in 
westward  group  velocity  due  to  the  eastward  jet 
(I  la);  the  opposite  is  true  for  the  westward  jet  ( I  lb) 
as  shown  in  Fig.  8b.  These  Doppler  shifts  are  sim¬ 
ilar  to  those  calculated  for  meridionally  uniform 
flows  of  the  same  intensity  (dashed  curves)  but  are 
much  less  pronounced  because  of  the  narrowness 
of  (I  la)  and  ( I  lb).  In  addition  to  these  Doppler  shifts 
there  is  also  a  small  frequency  offset  which  is  toward 
higher  values  in  the  eastwaid  and  lower  values  in 
the  westward  flows  (Table  2).  Frequency  shifts 


caused  by  the  antisymmetric  jet  (lie)  are  not  evi¬ 
dent  (Fig.  8c)  because  the  net  transport  across  any 
meridian  is  identically  zero. 

In  the  Pacific  profile  which  has  both  symmetric 
(Equatorial  Undercurrent)  and  antisymmetric  (South 
Equatorial  Cunenl- North  Equatorial  Countercur¬ 
rent)  components  of  the  mean  How,  Doppler  shifts 
are  virtually  absent.  This  is  because  the  symmetric 
component  is  weaker  than  the  eastward  jet  of  1 1  la) 
by  a  factor  of  2  and  because  the  antisymmetric  com¬ 
ponent  is  inefficient  at  producing  net  advection  of 
the  wave  fields  as  was  seen  for  the  jet  (1  Ic).  There 
is  however  a  small  frequency  offset  in  Fig.  fid  sim¬ 
ilar  to  that  observed  for  the  eastward  jet.  Note  that 
the  point  of  zero  group  velocity  font  •  I  and  2 
modes  in  the  Pacific  flow  field,  as  well  as  in  the 
more  idealized  fields  of  Fig.  2,  is  always  near  zero 
wavenumber. 

4.  The  Kelvin  wave 
a.  Method  of  solution 

A  perturbation  scheme  is  developed  that  exploits 
the  smallness  of  R  and  the  simple  expressions  for 
the  Kelvin  wave  in  an  ocean  at  rest  given  in  (8). 
We  first  expand 

•  ■tt7l«,t  R:tt,  +  . . . ,  (14a) 

v  -  o,  +  R«,  +  R*»,  +  . . . ,  (14b) 

p  m  p%  +  Rp,  +  RVt  +  -  -  .  ,  (14c) 

<r  •  <r,  ♦  Rtr,  ♦  RV,  +  . .  %  .  (I4d) 

In  addition  we  expand  A*  for  reasons  noted  in 
Section  2. 

A*  -  A,*  +  RA,»  +  R*Af*  +  . . . ;  V-f  fiD,.  (I4e) 

When  substituted  into  (I),  one  finds  to  zeroth  order 
in  the  Rossby  number  the  Kelvin  wave  of  (8). 
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VoLUMS* 


Fig.  I.  Dispuruou  curve*  for  the  mode*  n  •  0.  I.  2.  el  tmall 
wavenumber  and  m  the  pretence  of  tbc  eastward  jet  (a)  (he  west¬ 
ward  jet  (b>  the  antisymmetric  jcl  ic)  and  the  Pacific  profile  id). 
Solid  fine*  give  dispersion  curve*  for  no  background  flow .  pkmed 
symbols  define  curve*  with  background  flow.  In  lai  and  lb)  addi¬ 
tional  dashed  curves  show  dispersion  relation  w  hich  would  result 
Rom  lutitudinatty  uniform  background  flow  with  velocity  equal 
to  jet  maximum  velocity. 

viz-,  *»  “  Po  »  4b,  tr»  *  k.  Since  v,  »  0  the  first- 
order  equations  are 

-to** i  -  ye,  +  ikp ,  »  -IkUu,  +  »>,«„  (IJa) 
+  pa,  +  p„  -  0,  (15b) 

-i<Tt?  i  ♦  iku ,  +  B„ 

mi°'  +  75D.-kL')p' (,5c) 


The  right  side  of  (IS)  can  be  written  in  terms  of 
Hermite  functions  by  projecting  Uu,  -  Up,  m  U*, 
and  D,u,  »  D,*,  onto  the  set  (*.}: 

f>« *t  m  D,  £  d.t It,, 

d,  -  D."  J*  D.iMudy.  ()<•) 

V*.  -  If.*.,  f.  -  f*  U*,*.dy.  (16b) 
a*0  I  — 

Since  each  t|r„  on  the  right  side  of  ( 13)  will  excite  a 
different  response,  it  is  convenient  to  further  de¬ 
compose  the  CXR)  problem  into  the  series 


“l  *  *1#  +  2  «u. 

(17a) 

• 

»i  -  »»  +  2  v„. 

■-I 

(17b) 

m 

Pi  *  Pn  +  2  Pi.- 

*>l 

(17c) 

The  first  term  in  these  expressions  deserves  special 
attention  because  dio  forcing  will  drive  a  resonant 
first-order  Kelvin  wave  with  i'„  -  0,i/,o  »  p„  «  <ha. 
In  order  to  remove  this  spurious  resonance,  one 
must  insist  that  no  di„  terms  appear  in  the  forcing. 
This  is  accomplished  by  imposing  the  constraints 
<r,  ■ fjt ,  U8a) 

Ai*  •  rf«g7>.  (18b) 

The  first  of  these  results  is  a  Doppler  shift 
proportional  to  the  wavenumber  k.  the  Rossby 
number  R.  and  the  width  and  direction  of  the  back¬ 
ground  Row  through  }„.  The  second  adjusts  A.  the 
eigenvalue  of  the  vertical  problem  underlying  our 
two-layer  limiting  case.  The  part  of  the  first-order 
solution  proportional  to  rfi,  (zeroth  order  Kelvin 
wave)  is  now  of  arbitrary  amplitude  which  can 
be  set  to  zero  without  loss  of  generality.  For  other 
n  *  I,  the  equations 

-rov<„  -  yt>„  +  ikp„  -  -ikf, *„  (19a) 

-itr,v„  +  yir„  +  p,„  -  0.  (19b) 

-i<r»Pi.  +  iku i,  +  -  -ik(f.  +  J,)*..  (19c) 

can  be  solved  by  eliminating  u„  and  p„  in  favor 
of  p„  to  get 

fi.w  +  (»•’  -  **  -  kta,  -  >•*)», » 

-  -2 ikiS,  ♦  Mrf.X*.,  +  y*.).  (20) 
Making  use  of  the  identities  a,  «  k  and 

-  0.  (21a) 

,:,w 
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/  H  +  1  \*a  /  fl  \,n 

>*.  “  ^ — j— j  *.*i  +  (21c) 

it  is  readily  determined  that,  for  n  »  1. 

»w  -  «*.**-..  c.  -  IWH/,  +  tW.>.  (22a) 


The  usual  method  of  obtaining  expressions  for  u 
and  p  from  v  by  combining  continuity  with  the 
zonal  momentum  equation  fails  here  because  <r, 
»  k.  But  from  (19b)  and  (19c)  one  finds 

*»  “  +  r.*,  +  (22b) 

Pi>  *  <7«9.*«  +  (r.  +  Vid.)*.  -  (22c) 

where 


«• 


(/.  +  Hd,)(?..-ir-^,,j  -  !/w.. 


(23a) 

(23b) 


t.  “  -  W.  +  VW„)(-!yj-L),'\  (23c) 

Thus  to  O(R’)  we  have 


"  -  +  R  \  lq.-i  +  r.+s, (24a) 


Fio.  9.  Dijpcmon  curves  for  the  Kelvin  wave  in  the  presence 
at  various  mean  Dows  specified  hyiR.di.  Solid  line  <R  *  0)  cor 
responds  to  the  Kelvin  wave  in  the  absence  of  mean  currents. 


lions  for  the  opposite  sign  since  reversing  the  direc¬ 
tion  of  the  mean  flow  simply  changes  the  sign  of 
all  the  first-order  fields.  Thus,  only  one  set  of  cal¬ 
culations  needs  to  be  done  for  eastward  and  west¬ 
ward  flows  of  the  same  width  and  intensity. 


(24b) 


I)  Dispersion 


P  -  +  R 

X  f  fo.-t  +  r,  +  V4i/„  -  (24c) 

<r  “  <r«(  1  +  R/i).  (24d) 

A*  -  g’/)(D.  +  Rtf.),  (24e) 

with  q-i  -  q»  “  0.  In  order  to  standardize  results 
for  various  cases  we  will  fix  A*  for  a  given  c  as 
discussed  in  Section  2.  This  requires  simply  that  the 
zeroth-ordcr  mixed-layer  depth  be  defined  by 

D,  -  1  -  Rtf, 

which  is  a  function  of  U  through  d%. 
b.  Results 

Solutions  to  (24)  were  computed  in  the  presence 
of  mean  jets  of  the  form 

Rl/  -  sRf-*,  a  >  0,  (25a) 

for  which  the  coefficients  (/.,</,)  readily  can  be  cal¬ 
culated  as  shown  in  the  Appendix.  The  upper  layer 
thickness  associated  with  this  flow  is 

D  -  1  ±  (/Mo*-*.  (25b) 

Solutions  to  (24)  for  one  sign  of  (25)  determine  soiu- 


From  (16b),  (24b)  and  (25a).  the  frequency  shift  is 

R«t,  -  R/0A  »  ±Rk(a  +  I)-'*.  (26) 

The  magnitude  of  this  Doppler  shift  is  greatest  in 
broad  (small  a)  mean  flows  and  as  expected,  higher 
(lower)  frequencies  result  from  eastward  i  westward) 
flows.  In  the  inertio-gravity  range  (Fig.  8)  similar 
tendencies  exist  but  are  more  subdued  because  the 
phase  speeds  of  those  waves  are  much  larger  rela¬ 
tive  to  the  mean  currents. 

The  limit  a  «  0  (uniform  flow)  is  instructive,  if 
we  relax  the  geostrophic  constraint  (2b)  and  simply 
set  D,  -  I.  In  this  case  we  obtain 

»i  ■ 

so  that  the  total  solution  is  just  the  Zeroth-order 
Kelvin  wave  (8)  with  frequency  <r«(  I  ±  R).  All  terms 
with  summations  overn  is  (24a)-(24ct  vanish  in  this 
(somewhat  artificial)  limit,  as  can  be  seen  from  defi¬ 
nitions  (16a)  and  (16b). 

The  dispersion  of  free  Kelvin  waves  in  a  variety 
of  mean  flows  is  shown  in  Fig.  9.  To  the  approxi¬ 
mation  (24)  all  the  waves  remain  nondispersivc  with 
frequency  depending  on  the  width  and  intensity  of 
the  mean  flow  as  well  as  the  Kelvin  wavenumber. 
Moreover,  as  is  clear  from  (26),  to  each  of  these 
curves  there  corresponds  a  family  of  overlapping 
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curves,  e.g.,  the  curve  (n  «  3,  R  »  0.2)  corresponds 
to  curves  for  (a  «  8,  R  «  0.3),  (a  »  0.  R  »  0.1), 
etc.  It  is  noteworthy  that  the  dispersion  plot  for 
any  antisymmetric  Jet  coincides  with  the  line  R  -  0 
since/*  ■  0  by  ( 16b):  as  in  the  case  of  inertio-gravity 
waves  discussed  above,  antisymmetric  flows  do  not 
Doppler  shift  the  wave  fields.  Observe  also  that  this 
figure  for  Rl/  >  0  is  similar  to  Fig.  I  of  Philander 
(1979)  which  shows  the  numerically  determined  dis¬ 
persion  of  Kelvin  waves  in  similar  eastward  flows. 

2)  Meridional  structures 

We  truncated  the  summations  (24)  at  n  «  40  (20 
nonvanishing  terms,  due  to  symmetry);  this  is  a  very 
accurate  representation  except  for  exceedingly  nar¬ 
row  (a  *>  5)  jets.  Figs.  10- 12  show  the  resultant 
modifications  to  Kelvin  waves  of  various  frequen¬ 
cies  in  eastward  and  westward  Jets  of  various  widths. 
In  Fig.  10.  a  long  wave  interacts  with  a  narrow  jet. 
The  first-order  fields  u ,.  i>„  and  />,  take  the  form 
of  low-frequcncy  Rossby  waves,  i.e..  «,  and  p,  are 
nearly  in  geostrophic  balance.  The  total  fields  in 
Figs.  10b  and  10c  are  nearly  the  same  as  their  R  «  0 
counterparts  except  for  the  zonal  velocity  compo¬ 
nent  near  the  equator.  The  frequency  is  10%  higher 
(lower)  than  <r,  ■  k  in  the  eastward  (westward)  jet. 

Broader  mean  flow  will  have  even  less  effect  on 


the  meridional  structures  of  the  Kelvin  wave  in  this 
region  of  dispersion  space.  As  shown  in  Fig.  1 1  for 
jets  nine  times  wider  than  those  of  Fig.  10,  the  first- 
order  fields  are  much  weaker  and  the  total  fields 
are  virtually  indistinguishable  from  the  zeroth-oider 
curves.  On  the  other  hand,  the  frequencies  are 
shifted  more  for  this  mean  flow  in  accordance  with 
(24d). 

For  frequencies  that  fall  within  the  inertio-gravity 
range  a  >  1 .  the  first-order  corrections  change  char¬ 
acter  dramatically.  First,  from  (22a)  the  meridional 
velocity  is  a  linear  function  of  wavenumber  and 
hence  frequency  so  that  this  component  is  much 
larger  than  at  lower  frequency.  Second,  the  correc¬ 
tions  forced  by  the  zeroth-order  solutions  are  dis¬ 
tinctly  gravity-like  because  of  the  nearness  of  those 
waves  to  the  Kelvin  dispersion  curve  at  higher  fre¬ 
quency  (Fig.  I).  Both  these  effects  are  seen  in  Fig. 
12:  v t  is  about  10  times  larger  than  that  in  Fig.  10 
for  a  low-frequency  wave  and  the  components  of 
the  first-order  field  look  very  much  like  those  asso¬ 
ciated  with  an  n  »  1  gravity  mode.  When  combined 
with  the  zeroth-order  solution,  these  forced  correc¬ 
tions  produce  a  much  different  picture  of  the  Kelvin 
wave  in  both  eastward  and  westward  shear  flow. 
We  note  that  Philander' s  (1979.  Fig.  I)  numerical 
result  for  an  800  km.  3.3-day  wave  in  narrow  east¬ 
ward  shear  flow  (R  —  0.28)  is  very  similar  to  that 
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Pm.  10.  Meridional  structure  of  a  tow-frcquency.  hug  wavelength  Kelvin  wave  in  narrow, 
shear  Dow  (dotted  curves  in  (h)  and  id)  First-order  corrections  art  shown  in  lal.  Solid 
carves  in  fb)  and  (cl  are  the  local  ReWs  in  the  presence  of  the  indicated  mean  How;  dashed 
carves  are  the  Kelvin  wave  it  the  absence  of  mean  Row. 
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•  observed  in  Fig.  12b  which  is  for  a  900  km,  3.3- 
day  wave. 

S.  Discussion 

The  foregoing  analysis  has  shown  that  meridional 
velocity  fluctuations  associated  with  neutrally  stable 
incrtio-gravity  waves  advect  the  mean  currents  so 
as  to  produce  intermediate-frequency  meanders. 
The  strongest  of  these,  associated  with  the  mixed 
Rossby-gravity  wave,  are  similar  to  those  observed 
during  Phase  II  of  GATE.  Likewise,  Titov  11977) 
reports  on  meanders  of  the  Cromwell  Current  be¬ 
tween  2.5*N  and  2.3“S  along  166. VE  in  the  Pacific. 
Over  a  two- week  period  in  March  1972.  the  maxi¬ 
mum  speed  of  the  mean  current  migrated  about  the 
equator  with  a  period  of  0(  10  days),  lagging  (leading) 
the  meridional  velocity  by  about  a  quarter  cycle 
north  (south)  of  the  equator.  These  phase  relation¬ 
ships  and  symmetry  characteristics  suggest  a  mixed 
Rossby-gravity  wave  interpretation  like  that  used 
to  describe  the  GATE  meanders.  It  is  noteworthy 
that  these  two  sets  of  current  measurements  (the 
most  comprehensive  to  date  at  low  latitudes)  dis¬ 
play  the  same  type  of  meandering  behavipr.  for  it 
suggests  that  such  meanders  are  commonplace. 
However,  the  coincidence  may  be  fortuitous.  In¬ 
deed,  the  mixed  Rossby-gravity  wave  is  the  gravest 
antisymmetric  equatorial  mode  and  should  be 
strongly  excited  by  large-scale  meridional  or  anti¬ 
symmetric  zonai  winds.  Unlike  the  higher  order  in- 
ertio-gravity  waves,  though,  there  is  no  preferred 
response  to  wide  brnd  forcing  at  intermediate  fre¬ 
quencies.  implying  that  winds  with  highly  structured 
spectra  must  have  generated  the  observed  mean¬ 
ders.  Since  such  spectra  ate  atypical,  one  might  al¬ 
ternately  conclude  that  these  data  were  collected 
during  periods  of  anomalous  wind  forcing. 

Peaks  in  velocity  spectra  associated  with  incrtio- 
gravity  waves  other  than  the  mixed  Rossby-gravity 
wave  will  be  sharper  in  the  hemisphere  of  the  shal¬ 
lowest  mixed  layer  when  the  depth  increases  sig¬ 
nificantly  between  the  turning  latitudes.  Conversely, 
irrespective  of  the  pycnocline  topography,  sea  level 
spectra  will  be  little  different  from  those  predicted 
by  theories  without  background  flow.  This  may  have 
facilitated  the  work  of  Wunsch  and  Gill  (I9',6)  in 
modeling  4-day  and  5-day  sea  level  variability  as  the 
surface  manifestation  cf  baroclinic  inenio-gravity 
waves;  their  theory  included  no  mean  currents,  yet 
the  fit  of  observed  to  predicted  sea  level  amplitude 
distributions  was  good. 

Because  the  time  scale  T  of  (3)  is  proportional 
to  the  fourth  root  of  the  stratification  parameters 
(«,/>),  a  moderately  accurate  determination  of  the 
mean  density  field  fixes  T  and  hence  the  frequency 
rather  precisely  for  vertically  standing  modes  at  a 


given  zonal  wavenumber.  It  should  then  be  possible 
to  ascribe  significance  to  relatively  small  departures 
of  observed  frequencies  from  R>0  values  and  to 
interpret  such  shifts  in  terms  of  mean  currents.  For 
example,  Wunsch  and  Gill  ( 1976)  found  in  their  Pa¬ 
cific  sea  level  study  that  inertio-gravity  wave  fre¬ 
quencies  predicted  using  an  equatorial  /3-plane 
model  without  background  flow  led  to  frequencies 
that  were  uniformly  too  low  by  about  10%  compared 
to  the  observations.  Although  there  is  a  large  un¬ 
certainty  in  their  calculations  ( time  scales  were  com¬ 
puted  using  a  single  density  profile),  the  observed 
discrepancy  is  consistent  with  the  3-7%  higher  fre¬ 
quencies  expected  for  inertio-gravity  waves  in  the 
Pacific  profile  (Fig.  8d). 

The  low-frequency  Kelvin  wave  solutions  of  Sec¬ 
tion  4  are  very  similar  structurally  to  Kelvin  waves 
in  an  ocean  at  rest  except  for  their  higher  or  lower 
phase/group  velocities.  This  change  in  propagation 
speed  can  be  large  in  strong  unidirectional  currents 
and  thus  may  affect  significantly  the  baroclinic  ad¬ 
justment  times  in  the  tropics.  With  increasing  fre¬ 
quency.  the  meridional  velocity  of  the  Kelvin  waves 
increases  in  the  presence  of  mean  flow  and  they 
become  more  gravity-likc.  If  not  directly  wind 
forced,  these  waves  could  be  excited  by  gravity 
modes  in  western  boundary  regions  and  propagate 
into  the  open  ocean  as  varicose  meanders  of  the 
mean  currents.  Their  energy  contribution  to  sea 
level  and  velocity  fluctuations  may  be  comparable 
to  that  of  gravity  waves  in  the  same  frequency  band, 
but  separating  these  two  types  of  motion  requires 
wavenumber  spectra  as  well  as  frequency  spectra; 
at  present  such  data  are  lacking. 

The  basic  simplification  of  the  foregoing  analysis 
is  the  two-layer  model  with  its  single  baroclinic 
mode.  Unfortunately,  there  are  large  regions  of  dis¬ 
persion  space  where  a  single  vertical  mode  poorly 
describes  the  field  of  motion  (Philander.  I97gbi.  For 
realistic  stratification,  this  occurs  when  the  forcing 
is  confined  to  a  thin  surface  layer  and  when  its  fre¬ 
quency  and  wavenumber  arc  such  that  a  large  num¬ 
ber  of  vertical  modes  are  excited  nonresonantly. 
Under  these  conditions,  the  ocean's  response  is  bet¬ 
ter  described  as  a  sum  of  vertically  propagating  hori¬ 
zontal  modes  as  demonstrated  in  Wunsch  ( 1977b) 
for  periodic  monsoon  forcing.  It  is  important  to 
know  how  these  vertically  propagating  modes  inter¬ 
act  with  mean  currents  near  the  surface,  a  problem 
that  is  the  natural  complement  to  the  one  considered 
in  this  paper. 
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